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Heat Transfer and Reynolds’ Analogy 
in a Turbulent Flow with Heat Release’) 
By Raut R. Hunziker, Washington, D. C., USA.2) 


Introduction 


Pat has found polynomial solutions of REYNOLDs’ turbulent flow equations 
- expressing the time average axial velocity and correlation distributions in a 
circular pipe [1]?). OSBORNE REYNOLDs was the first who inferred a dependency 
between the transference of heat and momentum in turbulent flow [2]. Accor- 
ding to the ‘momentum transfer theory’ [2,3] these distributions allow a 
consistant semiempirical determination of the eddy diffusivity for heat 
ey = % &y, (4, 5, 6}. 

-The Reynolds analogy for a circular pipe can be formulated departing from 
the axial Reynolds equation and the temperature distribution equation [3, 7]. 
It follows that the nondimensional temperature and velocity profile are iden- 
tical if the fluid contains heat sources of intensity proportional to the constant 
axial component of the gradient of pressures and if the effective Prandtl 

number Pr* = Pr « of the turbulent fluid is unity. Extensions of the analogy 

for Pr* + 1 have been developed by PrAnprtr [3] and TAyLor [8] assuming 
a laminar layer close to the wall in which heat is transferred by conduction 
alone. 

This separation of the field into laminar and turbulent regions was further 
improved by von KArMAN, who introduced a transition of buffer layer in 
which both molecular and turbulent heat transactions were assumed [9, 10]. 
MarTINELLI has introduced the effect of molecular diffusivity in the turbulent 
core which is of special significance for the case of liquid metals [11]. 

Similar analysis allowing the representation of the radial distribution of 
temperatures at great distance from the pipe inlet have been developed by 
Lyon [12], SEBAN and SHIMAZAKI [13]. . 

The calculations based on the extensions of the analogy are usually com- 
pared however with measurements of heat transfer from fluids in which no 


sources of heat are present. 


1) The results presented in this paper were obtained in the course of research sponsored by the 
Office of Scientific Research of the ARDC, US Air Force under Contract AF 18 (603)—104. 

2) Reed Research Inc., Chief Mathematics Branch, Research Division. 
Present address: Convair, A Division of General Dynamics Corporation, San Diego, California. 

3) Numbers in brackets refer to References, page 314. 
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The solution of the boundary value problem for the temperature ction 
can be constructed also for the case of a distribution of heat sources within the 

fluid. It will be shown that the Reynolds analogy can be extended to a ee 
releasing turbulent flow, allowing an estimate to the approximation given by — 
the series solution for the temperature distribution. 


oC ent ve th hy 


1. Formulation of Equations q 


For the case of fully developed turbulent flow in a pipe of radius a defined by 
OSn=—S1, 050522, -wSE==S+00 


the average velocity components and fluctuations are assumed 
“a=0,-v=0, @=0), ¢=4) v =3, o-oo ee (1)%' 


along the cylindrical coordinates 7, 0, € respectively. - 

The absolute temperature at any point has a time average T anda fluctua-_ 
tion 7” of zero mean value. The mean value of the eddy heat transfer fluctuation 
ocu’' T’ is assumed to be a function of 7 only. According to the ‘momentum 
transfer theory’ [2, 4, 5, 6] 


1 oT 


=u T'(n) = at bas (2) . 
where €, = « €y is the eddy diffusivity for heat and 
a aa) 
Su) = aR i\lae ” 


defines the eddy diffusivity for momentum. The ratio « will be assumed 
a = 1 (4, 5, 6], vis the kinematical coefficient of viscosity and f,(), fo(7) are the 
Pai polynomials expressing the time average axial velocity distribution 
w(n) = w(0) f,(y) and correlation of velocities 


wu! wo’ (n) = (") w(0) fo(%) 


respectively. With the usual designation Pe is the Péclet number Pe = Pr Re, 
defined as product of the Prandtl number Pr = wc/k and Reynolds number 
Re = 2aw,,/v, where c is the specific heat, k the thermal conductivity and 
W» = w(0) the maximum of w(7) at the pipe axis A= 0. 

It will be considered that the fluid and wall of the pipe are maintained 
at a constant temperature 70, 7) up to a certain cross section € = 0 at which 
the temperature of the fluid at the wall is changed to a different temperature | 


— 
« eee? 
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TE1)=T, = const for the remainder (0 << +00) of the pipe. The 


distribution 7(&,7) is reduced to nondimensional form by the change of 


- variable 


U(E, 9) = Ee, ( 
se 
where 7) is a mean of the inlet temperatures T (0, n). 

It is assumed that the volume heat sources do not cause a convective effect 
of perturbation on the observed fully developed turbulent flow. The heat source 
distribution o is assumed a function of (€, 7) and the dissipation of energy due 
to viscosity and the variations of conductivity k, specific heat c, kinematic 
viscosity v and density o due to temperature will be also neglected. By analyzing 
the local energy balance and the corresponding boundary conditions the 


average temperature distribution 7(&,7) is defined by the solution of the 
following parabolic nonhomogeneous boundary value problem | 


0 dU dU . 

ay Ao) 3] — BO) ae = GEM). E> 0, OSnsT, | 

UE =0, 22 G0) =0, Ul) =lmUE»), | 
where 


Pe 7 o(&, 7) a® 
Ati) =n —2Prew)), Bey=—Z ah), Géa=— ee. | 


G6) = af, ln) =14 Ck + Can, fe) = Con — 28). | 


Since the axial conduction is assumed small as compared to the large axial 
convection, the term 7 (0?U/0&) on the left side of the first of equations (5) 
has been neglected. 

From the experimental data of LAUFER [14] for Re = 50-000 it results [4,6] 


C, = =0-32742, Cy= —0-67258, C,= 43-0445, n= 32. 


2. Solution of the Boundary Value Problem 


The application of the Laplace transform [15] leads to a function «(y, p) 
defined in 0 <j <1 for real positive values of the parameter p = —A, which 
is a solution of a nonhomogeneous Sturm-Liouville boundary value problem 


L{o(7)] = oy [A (m) w'(y)] — P B(y) o(n) = —B(n) Ul) + Ln, A) » (7) 
a’ (0) = () ’ o(1) = 0 , | 
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where I(y, ?) = 2{G(é, n)} is the Laplace transform of G(é, 7). The solution of 
equations (5) is found by applying the inverse transformation {~? and the con- 
volution theorem to the solution of equation (7) 


UE, 1) = 3 e-™* com() [bm + 5% (1 — em 4)], (8) 
=0 m: 


m 


having assumed G(é, 7) = G(y) and where /,,, and @;,(7) are the eigenvalues and _ 
eigenfunctions of the selfadjoint operator L[{w/(n)] of equation (7). 
The coefficients are defined by 


dm = | Bln) Uo(u) m()) 4M » m= [e0)0 m(1) 47 « (9) 


Hence we find that the heat flux distribution is given by 


ee Sein ) Soult eval es eel ae (10) 


or 
m= 


The system of eigenfunctions {o m(1) is complete in the sense of approximation 
in the mean and is orthogonal relative to the weight function B(y) => 0 over 
[0-1]. Since the eigenvalue problem is definite positive, R1tz’s method provides 
sequences {A}, {o”(n)} convergent to the eigenvalues 2, and eigenfunctions 
«;(y) as has been proved [16, 17]. The system of ‘coordinate functions’ 


twin} = {1 — 9}, @ =0,1, 2,...), 


is complete and linearly independent and satisfies the boundary conditions. 
Thus the set 


Wm) =S) as yili), ag Uy eee 


is admissible. The extremal condition of the functional of which L{w(y)] = 0 is 
its Euler equation is given by the linear system 


Dd, % (li; — 4M, 3) = 0 CF Roe Ui eae 


j=0 
where 
1 
(2 + 4) (2 + 9) 4 f Gq (l= ht ge eer 
[==> —— —q Pry (2 +7) (2 i: 3\ | 
y 447+] Be) ee 2C,+2n€, 72%) dy, 
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, ay iRe fl 1 1 1 
eter: 5 Nes eee Sag, easy 4474 : 
| 1 1 1 1 
| c.| Ea Se 
- pata Say “Os ee 
1 t 1 
C.| pats 
ve SES Pee | 2n+4414 


1 
eee 

The ‘secular’ equation of degree m, P(A) = det|I,; — 2.M,,;| = 0, has m real 

positive solutions A (k = 0,1, 2,...,m), which are upper bounds of the 

corresponding eigenvalues. For each one of the Am”) each homogeneous system 
~ has a solution {a)”)} which is finally characterized by the normality condition. 

Assuming a constant distribution o(&, 7) = o, the coefficients in equations (9) 
are reduced to 


Pe 1 1 A {oe 1 1 
“Sra cit eee A See Ps 
Om = 2 a4 E tape Oly saa) + lsaaa LEP REAIL 
1 1 
(m) 
a Se 
(ee ay a2 a : (; #4) 
be : h (Ty — Ty) 


The temperature and local heat transfer distributions are obtained along the 
whole extension of the wall and not only in the radial direction of a long pipe 
as in other solutions [11, 12, 13, 18]. 


3. Extension of Reynolds Analogy 


OSBORNE REYNOLDS [2] was the first who inferred a dependency between 
the transference of heat and momentum in turbulent flow. As is known, 
L. Pranptt [3] applied REYNoLps’ idea of proportionality between the velo- 
city and temperature field near a heat exchanging surface by comparing the 
differential equations of the velocity and temperature fields. In this case of 
turbulent flow in a pipe the Reynolds differential equation along the axis is 
(see equations (19) and (42) of reference [1)) 


d } em(n) ]\ _ oz : 
a v £,() [1 at eal eat Oy,” (12) 


where 0/02 is the constant derivative of the mean pressure p. By equation (8) 
it is assured that 
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hence the solution given by equation (8) of the boundary value problem 
[equations (5)] for the temperature satisfies asymptotically the ordinary dif- 
ferential equation 


aE ult | 
5 SR Nae 
d he v bores (7a 13 
1 Oe a ee ee i 
assuming 
o=06,4+6,; y= (14) 


Since the boundary conditions of equation (12) and (13) are identical, 


f(co, 1) = U(co, 1) =0, f’(cv, 0) = U’(o, 0) = 0, 


it follows 
re Ath» aD 
if and only if 
PA =o Pre LS oe 6, =i, == COnsh, 
op 
Fe on = ae (16) 
kl, Tie eee 
Wy h (TF T.) op 
‘ 0 fag Oz 


which leads to the following 


Theorem: Consider a fully developed turbulent flow containing volume heat 
sources of constant density o, = 0, + 0, as defined by equations (14) and (16) with 
constant temperature T,, at the walls of the conduit and with mean velocity given 
by w(n) = Wm f(y). The asymptotic distribution of nondimensional temperature in 
the fluid ts 

id eal SS 5 SINS EW ertine WPS eR eS et (17) 
TT 


w 


ay and only. tf Pr*® = «% Pr = 1. 


When o, = 0, that is when y = 0, the theorem reduces to the original form 
of PRANDTL. 


4. Conclusions 


This generalization of ‘PRANDTL’s heat source theorem’ permits us to 
obtain an exact asymptotic distribution of temperatures. If the asymptotic 
nondimensional temperature profile U(7) for Pr* = « Pry = 1 is approximated 


“ty 


% 
=a 
3 
> 
i 
\ 
, 
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by the series solution [equations (8) and (9)] taking into account the presence 

of virtual sources [equations (16)] within the fluid, it can be compared with 
the exact solution [equation (17)}. 

Hence it is possible to estimate the error performed by using Rrtz’s method 

of approximation and the consequent series expansions for Py* = 1, z=oo. 

The same comparison and application holds in the case o, = 0 or » = 0. Since 


_ the shear stress at the wall t, satisfies the conditions 


= fee a op ae i tam! dw 
0 oR Tamed a Si pigs) 


Free hla nal 


oT = il y 
goo) =—h(-) = th (Th— Te) C, Re. (19) 


This equation expresses REYNOLDS’ analogy of heat transfer and fluid friction 


_ which is valid if and only if the fluid contains heat sources of density 


"O70, + 6, = 0, (Ley) and Py* = Pre =1. 


If the density of real sources o, is zero it follows that the density of sources 
for which REYNOLDs’ analogy applies is ¢, = 0, given by equation (16) as 
directly proportional to the pressure derivative along z. Similar results hold for 
turbulent flow between parallel plates [6]. All these formulae and consequences 
hold for the particular case of laminar flow, where /,(7) is represented by a 
parabola, and equations (5) reduce to the Graetz-Nusselt equations [19, 20]. 
With consideration of a rectangular velocity distribution, HARRISON [21] has 
calculated different solutions for the temperature with and without the axial 
conduction term 7 (0?U/0&). He concludes that the effect of neglecting the 
axial heat conduction term is not significant for Pe = 40 which will be pro- 
visionally considered as the lower limit of applicability of this analysis. This is 
the only limitation of the present analysis which may be applied to any incom- 
pressible fluid in fully developed turbulent flow, and in particular to the heat 
releasing turbulent flows which occur in nuclear reactors. The design of circu- 
lating nuclear fuel — reactor systems requires information on the distribution 
of temperatures in the fissioning fuel. The heat source distribution function in 
the moderator sorrounded channel of a liquid fuel reactor may be considered 
as G = G(n) in equations (8), (9) and (10). The values of 2 and a™, (=m 
=0,1,2.) have been calculated for five Prandtl numbers in the liquid metals 
region (0-002 < Pr < 0-030) at a Reynolds number Re = 50-000 and these 
values can be used for the evaluation of equations (8), (9) and (10) [4, 6). 
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Using equation (10) for G(j) = 0, « = 1 a closer agreement with experimen- | 
tal data [22, 23] is obtained than that indicated by the Seban-Shimazaki | 
equation [4, 6, 13]. 
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bi 
Zusammenfassung 


; Es wird die Warmeiibertragung von und zu einer voll entwickelten, warme- 
_ abgebenden, turbulenten Strémung in einem Rohr analysiert. Durch die Anwendung 

der Laplace-Transformation wird die partielle Randwertaufgabe der Temperatur- 
_verteilung in einem Rohr von kreisférmigem Querschnitt auf eine gewdhnliche 
_ Sturm-Liouvillesche Randwertaufgabe zuriickgefiihrt. Die entsprechenden Eigen- 
_ werte und Eigenfunktionen kénnen nun mittels des Ritzschen Verfahrens fiir ver- 
_ schiedene Werte der Prandtlschen und Reynoldsschen Zahlen berechnet werden. 
Die Verteilung der Temperatur und des Warmeflusses entlang der Wand erhalt 
man in Form einer Reihe, deren Glieder diese Eigenwerte und Eigenfunktionen 

enthalten. Durch Verallgemeinerung des Prandtlschen Warmequellen-Satzes wird 
_ die Reynoldssche Analogie der Warmeiibertragung und der fliissigen Reibung fiir 
_ eine warmeabgebende turbulente Stroémung erweitert. 

Die genaue asymptotische Warmefluss und Temperaturverteilung, die man 
erhalt, gestattet eine Abschatzung des Fehlers, den man bei einer Anwendung der 
Ritzschen Methode und der daraus folgenden. Reihenentwicklung erhalten wiirde. 
_ Diese Analyse kann ftir die warmeabgebenden turbulenten Str6mungen in Atom- 

reaktoren angewendet werden. 


Ww 


— 


(Received: April 15, 1958.) 


On the Phenomenon of Brake-Fading") 


By FREDERICK F. Linc and EDWARD SAIBEL, Troy, New York, USA.?) 


Introduction 


The importance of knowing the braking-time and braking-distance required 
to bring a vehicle to rest is apparent. Such calculations however, are quite 
difficult due to the variation of the coefficient of friction with temperature. 
This variation brings about a phenomenon known as brake-fading. 

In general, the temperature rise due to the friction causes a decrease in the 
coefficient of friction. The decrease in friction in turn affects the generation 
of heat and in this way an interaction takes place between the heat generated 
and the dynamics of the problem. The net effect under these conditions is to 
increase the time necessary to bring the system to rest. 


1) This work was sponsored by the Office of Ordnance Research, US. Army. 
2) Rensselaer Polytechnic Institute, Department of Mechanics. 
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It is the object of the present paper to indicate how, for a given relation — 
between coefficient of friction and temperature, this interaction may be © 
considered and the actual braking-time and braking-distance calculated. 


Theory 


Let V be the relative sliding velocity between two bodies which constitute 
a brake system, M be the effective mass to be decelerated, f be the coefficient 
of friction between the rubbing materials, A be the brake contact area and 
b be the normal pressure between the bodies in contact. The law of motion 
is simply 


dv 
Mia sig heehee (1) 


where R is drag or windage force which will be neglected in the present work. 
If M,f, A and # are constants, integration of equation (1) yields 


- E Ae; 
eee ae (2) 


where V, is the initial velocity. It follows that the braking-time, i.e., the time 
forV = 0,isV, Mij p A. 

In many practical situations however, / is known to decrease with the 
increase of interface temperature. This is caused by heat generated during the 
braking process. Thus the braking-time cannot be computed from equation (2). 
The phenomenon of the increase of braking-time due to the aforementioned 
temperature dependence is known as brake-fading. 

In what follows it is proposed to examine a particular type of brake-fading 
in which f is of the form 


f=f,—06), (3) 


where @ is the interface temperature above ambient. Also M, A and # are 
constants. Equation (1) can be expressed 


Vav,— 4? fpoey ar. (4) 


Equation (4) can only be used if 6 (t) is known. Now, consider the case in which 
the brake system consists of one insulating body and a conducting body. The 
more general problem involving conduction of heat in both bodies will depend 
for its solution on the knowledge of how the heat generated will be divided 


Evol. IXa, 1958 On the Phenomenon of Brake-Fading Sif 


a a 


_ between the two bodies. Up to now no rigorous solution of this problem exists. 

_ However an approximate solution for a quasi- stationary state has been given 

_by Brok [1]*). Since the object of this paper is principally to show how the 

~brake- -fading problem can be solved we refrain from treating this case because 
4 fof the extra complexities this method of partitioning entails. 

Let the area of contact be large so that a single interface temperature will 

_ describe the temperatures at all points over the surface. Thus the problem can 

| "be assumed a one-dimension flow problem in the direction perpendicular to 

the interface, z. For short duration of braking, the body can be considered 

semi-infinite. Mathematically, this situation is described by the heat equation 


+ 


076 1 00 


Oa oe Ob. (5) 


| where « is the thermal diffusivity. The initial condition is 6(z, 0) = 0. The 

boundary conditions are 6(0o, t) = 0 and 06/dz (0, t) = q, where q is the heat 
input. At the interface, z=0, g=/#V. It is known that the integral [2] 
4 satisfying equation (5) and the stated conditions is 


é — , 
ye Vt ry 


i i V 2/4alt—t’ dt’ 
ee ee | PV | 2/sa(t-2) (6) 


where y is the density of the conducting materials and c is the Specific heat. 
- But, of interest here is only the interface temperature 6(0, t) = G(t ). As z +0, 
_ equation (6) becomes 


ee aa ye 7 
10 =F é = | (7) 


Thus the brake problem is reduced to solving equations (3), (4), and (7) simul- 
taneously. 
Introducing the dimensionless groups: 


is Vo M bce’ t 
= ae Se = —, i= —, 
as Vo?’ ie op A’ $ to f fo 
ve? |/ MMP he mets 
= -~ —= — as, ==>, 
pa Pe : A= e p 


3) Numbers in brackets refer to References, page 321. 


a 
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equations (4), (3), and (7) become 


=i-fueya, pele1-aTe, Te—f MME 
o(é) =1 | me) ae, wlé) ADEs ene J ae (8) 
Defining 
' o(&’) 
Is , ae aa ie | 
Eth gee 
the last two of equations (8) yield 
we) =1—2f We) KEE, 8) ae", (9 
0 


which is a Volterra integral equation of the second kind. The solution can be 
expressed in the form of a Neumann series 


lee) 


ue) = bode OayS (10) 


where 


Semis) a) 164% (a swale, ow eta. iq wie w6 te' le) ee ae 


Note that the series is an alternating one and its convergence depends on 4 and 


the range of €. The order of magnitude of 4 in an actual problem is discussed 
below. | 


The first of equations (8), and (10) yield 


S 


v(é) =1— | h + 3 (ayn Fe) | ae! 
or : 
of) =1-£- Yaomey, (11) 


niin, 


ie aa ia 


Laue a 
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where 


See 
GY = Hi J Ke, &") de" ae’, 


Cree 
G® =f KE, 87 GME" ae ae, 
06 


OO OOP OC ee ee Ca het CMOMCME Tey 


ret 
GeS [ [Ke E") Ge-ny(eEn) dé" deel 
00 


om - 
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h 
Ti e Chr —> 
es wm Figure 1 


Plot of dimensionless velocity v versus dimensionless time &. 


~ Since 


=, o(€") 
i ee 

(5, 6") Ea? 
equation (11) is an integral equation in the unknown v. This can be evaluated 
by an iteration process. The number of iterations needed depends on the 


parameter “ > 
a Vel pe P fo 


f= 29 t mad 


For 4 = 0, equation (11) degenerates into (2). 
Estimating from automobile data on weight, dimensions, pressure, thermal 
and material properties, and using a decrease of 25% in coefficient of friction 
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> 
in a temperature rise of 700° F, leads to a value of A of the order of one tenth. 
Thus rapid convergence is assured. : : 


Define v,, as the result of one or more iterations considering terms uP to 
those with multiplier 2”. Thus 


eget 


. 
. 12) 
ai 7 ( 
V=0,+77 ie (1 - 35). | ; 
1024 32 64 600 a 672 2 
oS Se eee neg 2 €3 
g—= tA 2835 [1 55 é+ (a5 me) aes I. 


Figure 1 shows plots of v versus & for several selected values of A. Figure 2 
shows a plot of the braking-time versus*/; also the percentage increase of time 
required for braking over the ideal case f = f, = const or 2 = 0. For 4 < 1/2 
the percentage increase has already exceeded 50%. 


2:0 T 


= 
(=) 
Qa 


% increase of braking time over ideal case of f=f,= const.—=— 


— 
oO 
i=) 


| 
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oOo 


<= 
fs 0 
x 
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= 
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Plot of braking time and percent increase of braking time versus the parameter A. 


Stopping distance may also be calculated. By definition 
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rN ear aay aTee 


i Let us define a dimensionless distance wy such that 


_ tsd 
; oe Vaitg 
4 
é where S is the actual distance. Thus 
| 1. ds 1 ds 
4 v= = 
‘ ew dbase Vo ee 
. Or 
| dy 
Co we . (13) 


_ From equations (12) and (13) the dimensionless braking-distance Wn for the 
n-th iteration is 


; 
; 


i Yn =] 2, (6) a8, (14) 


_ where &,,, the dimensionless braking-time for the m-th iteration, can be obtained 
from Figure 2. 


Conclusion 


The foregoing method of approach is applicable even for more complicated 
relations between coefficient of friction and temperature. The calculations 
_ would probably be considerably more complex and recourse might have to be 
made to machine computation. If such facilities are available, a more realistic 
model of the heat transfer process might well be considered. This could involve 
both bodies in contact being conductive as mentioned above, the actual 
geometry of the bodies taken into consideration, consideration of heat loss 
(forced or natural) to the surroundings, amongst other possibilities. 
Should materials be used for which the coefficient of friction increases with 
temperature, the calculations for what may be called anti-brake-fading can 


also be accomplished. 
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Zusammenfassung 


Das Problem des Bremsschwindens wird fiir einen linear mit der Temperatur 
abnehmenden Reibungskoeffizienten untersucht. Die verwendete Methode lasst 
sich auch auf kompliziertere Zusammenhange zwischen der Temperatur und dem 


Reibungskoeffizienten anwenden. 


(Received: May 16, 1958.) 


Exact Transient Solution 
of the Buried Line Source Problem for an Asymmetric Source 


By Maninpra Mitra, Calcutta, India?) 


Introduction 


The problem of the buried symmetric line source in a semi-infinite elastic 
medium has been considered by NAKANO [5]?), Lapwoop [4] and GARVIN [3]. 
NAKANO considered both an oscillating and an impulsive dilatational source 
but Lapwoop and GARVIN considered only the latter case. LAPwoop applied 
Fourier integral methods and obtained approximate expressions for the solu- 
tions while GARVIN used CAGNIARD’S method to obtain the exact solution. In 
this paper, the latter method has been used to determine the surface displace- 
ment components in the case of an impulsive dilatational source buried in a_ 
semi-infinite medium. The impulsive character of the source is different from 
that considered by GARVIN. 


Problem 


The line source is assumed to be situated at (0, 2) in Figure 1, the homo- 
geneous, isotropic elastic medium lying in the region y => 0 (Figure 1) and the 


yy 


Figure 1 


Hlastic medium with position of source indicated. 


1) Jadavpur University, Department of Mathematics. 
2) Numbers in brackets refer to References, page 331. 
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_ medium to be at rest before ¢ = 0. The surface displacement components will 


be obtained for a source represented by the wave-function qo, given by 


9l7,9,4=0, t< — 


a cos 6 y2 Y Y 
es 2 
:> Y /: : eS Fol (1) 


a2 


a cos 6 ie Oe Y 
= £9 new F. Sram 


where (7, 0) are indicated in Figure 1, « is the compressional wave velocity in 


the medium and q, satisfies the wave equation 


op op 1 dp 
Sar Tiel ua (2) 


; everywhere except for r=0, and ¢=+7/a, ¢=i,+7/«. The radial and 
_cross-radial displacements (u,)9, (%#g)y due to the source alone are indicated in 


4 


_ Figures 2 and 3 respectively. Since the radial displacement consists of two 


sudden jerks in opposite directions, each followed by a gradual incomplete 


recovery, the impulsive character of the source is different from that con- 
sidered by Lapwoop and Garvin; in the case considered by them, the 
radial displacement consisted of a single sudden jerk followed by a gradual 
- incomplete recovery. 


Up, 
a cosOty 
pe 
rid Lr —>t 
Sf +h 
a la 
Figure 2 
Radial displacement due to the source alone. 
Ue, 


asin Ot 
pe 


2 
ato 


Q|> 


Figure 3 
Cross-radial displacement due to the source alone. 


ZAMP IXa/21 


324 MANINDRA MITRA ZAMP 


The components (w, v) of the displacement along the x- and y-axes (Figure 1) 
are given in terms of g and another wave function y by the relations 


1 


a Pe ie EP ign Pex Ts 3 
ua 2 ye (3) 
where 
o2y o7y 1 o7y 2 
ox i oy? “pe . or ( ) 


and f is the shear wave velocity in the medium. The stress components yy, x 
are given by 
—~ Ou 


ou — Ou ou 
yO) ee ES eee =n (= a x) 


(4) 


4, @ being the coefficient of rigidity and density respectively of the medium. 
The problem is to obtain a solution 


p(x, ¥, t) = olx, ¥, t) + pi(x, ¥, 2), p(x, y, t) = yylx, 9, 2), (5) 


where @ satisfies equation (2) and yp satisfies equation (2’) and which are such 
that (i) the stresses at the free surface y = 0 of the medium vanish and (ii) 
the displacement components (w, v) tend to zero as y > +00. 


Solution 


Defining the Laplace transform /(x, y, #) (for real positive values of #) of 
a function f(x, y, t) by the relation 


Hx, ¥, B) = fe? f (a, y, t) at 
0 
and denoting all transformed quantities by superimposed bars, we obtain from 
(2), (2’) the differential equations 
7 ge 
axt + aye Soa P (6) 


Oy ee ee 
ox ; oy? Bp yp . (6’) 


We now proceed first to obtain the Laplace transform @, of @) and then to 
obtain an integral representation for my for0 <y < h. 
It is easily verified that 


Yo= h(p) Ky (22) cos 6 , 
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_ where 
7 > —Ply 
fy(p) = fe? H (tt) dt =, 


as (t) being HEAVISIDE’s unit function. 
For @ satisfies equation (6) and K,(f r/«) is the Laplace transform of y(é) 
given by (ERDELY1) 


and so 


which is identical with qp. 
Now 7 cos@ = x, andr sin@ = h — y (Figure 1) 


and 
® .— (h—y) VES Ble ar : 2 
ae kink xdk=—_*Pl8 __ x [ P Via—ye + Al, 
J fk? + p3/o? V(A— vy)? + #2 o% 
: (ERDELYI) 
= a4 cos 6 Ky (=| ; 

we have 
ae fk oe pret Om 3 
pa= > Alf) / a2 sink dk, OS¥<h, (7) 


ern al Pk ple 
We assume that 


n= / BQ) er" sinkx dk, Pi= it C(k) e7'" cosk x dk, (8) 
0 ) 
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| 

where 7, = Vike + p/p? p?/p2. It is easily seen that 7, yy satisfy « equations (6), (6’) | 

respectively and that the contributions to w, v from @,, y, tend to zero as | 

Ver OO. ! | 
Using relations (5), (7), and (8), we obtain 


a= [|- a h{b) ty) np _  B(R) en? — 7, C(B) es] cosh x dh , 


Pn 
i ; (9) 
v= 12 oa ey + Bh) np ew?" + bC(R) e-? "| sink x dh 
0 


Using (4) and the conditions xy=yy=Oat y=O0s0 that xy= yy=0 
at y = 0, we obtain 


(6 + 2 RB) BYR) + (2 BP hy) CR) + 2" (p24 2 Be” f(A) = 0, 


(2. Bk ny) Blk) + (62 + 2B?) C(k) —- 222 ch F,(p) = 0 


P 
so that 
(0) Scar Sater agp 
CA ia OS, is See Bye pa ee ae 
Bi) = — Fe PAO) — Say PAE 


where A (k*) = (p? + 2 Rk? B?)? — 4 kh? Bt ny 7. 
Substituting these values of B(k), C (k) in (9), the following expressions are 
obtained for u,v when y = 0 


u(x, 0) = —4 pf, (p A | ja Be ep cosk x dk , 


U(x, 0) = 2% + he) {= a B? (p® + ey = Mis to Po) a ae 
0 


where vy is the y-component of the displacement at (x, 0) due to the source 
alone 


u(x,0)=—4p f(p) Re |/- ot cae? eh np—ikx a| 


0 


v(x, 0) 2 Up 4 ie Im / a B? RS (p? + mire — 2175p B?) eh Np—thx | : 
0 


” 
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ies hY1l+2+iz% 


Ewe eg” a 


i 
in the above expressions we get 


u(x,0) = Sees ae Re | 5 Bee e | 


0 


dy ee Tin [Be 2) onal, 
0 


~ where 
_— BY1+ 922 
A en 
1 222) _2y/1 + 2 232 


A(y, 2) 


pand A(y, z) = (1+ 2° 2’)? — 48 2V1 + 2V1 +92 22. 

4 The integrands have branch points at z= +7, +7/y, and simple (Ray- 
_ LEIGH) poles at z= --7 «/V,, where V, is the velocity of Rayleigh waves in the 
medium. As 

a 1 


Vr yas 


Ve<B, 
We shall now consider the transformation 


hylt2+izx 
= V ae : (10) 


4 i 


_ As z moves along the real axis from 0 to +o, ¢ moves along the curve AB to 
infinity where AB asymptotically approaches the straight line 


Oo 


z being a real parameter. As z decreases from + 7 through the \ value —ix|Vh2+ x2 
to — i, t increases along the real axis from — x/ax to Vn? + x?/a and then 


decreases through real values to x/x. For imaginary values of z outside the 
integral (— i, + 1), ¢ lies on the curves J,, J. The branch points and the 


Rayleigh poles are shown in Figure 4 (GARVIN). 
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mi cn er a aa a 


tplane 


Figure 4 
t-plane. * branch points. 
The integrals 
Cc 
dz 


Cc 
; 
ii Fi(v, 2) em dt and f Flv, 2) e-?t Fd 
B B 


tend to zero as BC moves to infinity. Since there are no singularities of the 
integrands inside the contour ABC DA (it should be noted that the portion 
AD of the contour is approached from the upper half-plane) 


foe) 


00 B 
; dz pee d 

i F,(v, 2) e-?* dz = / F,(v, 2).e7?# = JAS " Fi, z(t) enpt a dt 

0 A 


hia 


and similarly 


co : 
/ Fy(vya)er?! da = | Fv, 2(t)] e-?* — dt. 
7 hia 


Thus by inspection it is seen that 


Re x0: AN oe dg a 
Im | [ne ayes | 


and 
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: are the Laplace transforms of G,(¢) and G,(2) respectively, where 


EN =0, Oana, 


= Re [R(>, 2} 9 |). 2<:, 


(11) 


and 


is obtained by inverting (10). 
Now for 


_z is a pure imaginary such that |z| < 1, while dz(t)/dt is a pure imaginary. 
Hence F,(v, z) dz/dt is a pure imaginary and /,(y, 2) dz/d¢ is real for 


eA sachs 


We can, therefore, replace h/« in the inequalities in (11) by Vine + x2 /ot. 
Since 


p hip) = ert = | e-P8 (b— te) dt, 
0 
6(¢) being Drrac’s delta function, 
alr = AP fal 16 Aaa 
Baie =f 3 


=-*"G (-4), tyo<t. 
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Similarly, we obtain 
OT, Oe Age Peele: 


4 y 


From (1) and (3) we obtain 


Ge 
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; The equations (12) show that the horizontal displacement u(x, 0) at the 
surface begins a period ¢, after the beginning of the vertical surface displace- 
ment v(x, 0). 
I am grateful to Prof. B. B. SEN for his kind help and guidance. 
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Zusammenfassung 


Es wird die Verschiebung an der Grenzflache eines homogen-isotropen elasti- 
schen Halbraums, erzeugt durch die stossartige Dilatation einer eingetauchten 
Quellstrecke, exakt ermittelt. 


(Received: March 17, 1958.) 


Comportement d’un cable porteur-tracteur 


Par MARIE MARCEL NICOLAS, Paris, France?) 


On appelle cable porteur-tracteur tout cable pouvant se déplacer sur des 
appuis entre lesquels il est suspendu, tout en entrainant des charges qui lui sont 
accrochées et dont il supporte intégralement les poids. 

La détermination des formes diverses que peut prendre un cable porteur- 
tracteur et celle des tractions qui leur correspondent suscitent des problémes 
qui procédent tous plus ou moins des données schématiques suivantes: A et B 
sont les deux appuis successifs d’un cable pesant, sur lesquels il peut glisser 
sans frottement et agencés de telle sorte qu’en l’un d’eux, A, une traction 
constante est exercée sur le cable tandis qu’en l'autre, B, la traction du cable 
est appropriée 4 la longueur de cable suspendu ainsi qu’aux forces exercées 
entre A et B en des points fixes du cable, c’est-a-dire solidaires de celui-ci. 

La longueur et la forme du cable porteur-tracteur ainsi sollicité et la traction 
appropriée qu’il convient d’exercer sur lui au point B dépendent de la traction 
constante exercée en A, de la longueur et de l’inclinaison de la sustente A B, 


du poids du cable par unité de longueur, des poids des charges qui y sont fixés 
et des emplacements des points d’attache dans l’espace; elles dépendent encore 


1) Secrétariat d’Etat aux Travaux Publics, aux Transports et au Tourisme. 
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des trajectoires et des vitesses qu’on veut imposer a ces points d’attache, et par — 
conséquent du mode d’attache des charges. 

La confrontation de données aussi complexes ne peut étre utilement entre- 
prise tant que l’on n’a pas étudié le comportement d’un cable porteur-tracteur 
4 l’arrét, observation faite que cette éventualité peut toujours se présenter dans 
‘la pratique; observation faite encore que cette éventualité est aussi celle de 
toutes les installations (ponts suspendus par exemple) oti des cables fixes sup- 
portent des charges qui leur sont accrochées. © 

Telle est l’étude que nous allons exposer, en en simplifiant encore les données 
puisque nous ne traiterons que le cas d’une charge unique, appliquée en un 
point fixe, C, du cable. 

Plus précisément les questions auxquelles nous apporterons une réponse 
sont les suivantes: 

Etant donné deux appuis successifs, agencés comme nous I’avons dit, quels 
sont les points de l’espace ott une charge peut étre accrochée au cable? Quelle 
est la valeur de cette charge et aussi celle de la traction appropri¢ée en B en 
fonction du point de l’espace ot se trouve le point d’accrochage? A quelle 
condition peut-on disposer pour une charge donnée d’une ligne continue de 
points d’accrochage entre les deux appuis? Dans quelle zone de l’espace se 
trouvent ces lignes continues ? 

Dans nos raisonnements toutes les forces, et notamment les charges, sont 
exprimées en longueurs de cable dont les poids égalent les forces ainsi mesurées. 
Dans ce systéme de mesures l’ordonnée canonique [y = m cosh («/m)] d’un point 
du cable (dont la forme est celle d’un ou plusieurs arcs successifs de chainette). 
mesure donc la traction dans le cable en ce point. 

Enfin pour faciliter l’exposé, nous supposerons toujours que l’extrémité 
contrepesée A est figurée a gauche de l’extrémité B. 


1. Points de espace ot une charge peut étre accrochée 
au cable porteur-tracteur 


Sous l’action d’une charge unique le cable prend la forme de deux arcs de 


chainette, de deux caténes, AC et CB lesquelles ont méme module. S’il n’en 
était pas ainsi les deux tractions du cable, tangentes a ces caténes de part et 
d’autre du point C, ne pourraient pas avoir deux composantes horizontales 
égales et opposées et l’équilibre statique serait impossible. 

Si Aa est la traction constante exercée en A sur le cable, la caténe de sus- 
pension AC se trouve donc sur l’une des chainettes isosthénes du faisceau noué 


en A, et ayant comme base commune l’horizontale tracée A une distance A a 
en dessous du point A. 
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La chainette-support de la caténe de suspension AC renferme le point B 


_ dans sa zone concave. En effet, l'autre caténe de suspension CB, étant un arc 
appartenant a une chainette de méme module que celui de la chainette-support 


de AC a peut étre considérée comme ayant été découpée dans celle-ci suivant 


Varc C’B’ et avoir été transférée en CB par une translation définie par le 


vecteur B’B, ou son équipollent C’C (voir figure 1). 


(S) 


Figure 1 


Détermination des deux caténes de suspension afférentes 4 un point de charge C et de la charge 
pouvant y étre suspendue. 


Si B n’était pas dans la zone concave susdite, la translation C’C ferait se 
joindre en C, deux arcs de caténe ayant leurs tractions dirigées vers le bas, 
donc ne pouvant s’opposer qu’a une force verticale soulevante et non pas a 


une charge. 
L’appui B, devant étre dans la zone concave de la chainette-support de la 


caténe de suspension AC, il en résulte que celle-la ne peut étre l’une quelconque 
des chainettes isosthénes du faisceau A, et qu’elle se trouve nécessairement dans 


l’éventail compris entre la chainette-support de la caténe lége tendue, AB, et 


celle de la caténe lége lache, AB. 

Chacune des chainettes-supports est donc tangente a la courbe enveloppe (®) 
des chainettes isosthénes du faisceau noué en A et est complétement déterminée 
par son point de tangence?). Aussi désormais l’expression «chainette (~)» nous 
servira a désigner celle des chainettes-supports dont est le point caracté- 
ristique, c’est-a-dire le point de contact avec l’enveloppe (®). 


2) M. M. Nricoras, Cable de téléphérique supportant wne charge, Ann. Ponts Chaussées 1955, mai- 
juin. 
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Ces points g sont compris entre les points gj de la chainette-support de la 


caténe lége tendue AB (module: mo) et yo de la chainette-support de la caténe 


lége lache AB (module: mj) (voir figure 2). 
Parmi ces chainettes-supports il convient de distinguer celle dont le point 


caractéristique @ est a l’intersection de la sustente AB prolongée avec la courbe 


Figure 2 


Eventail des chainettes-supports susceptibles de servir Ala suspension d’une charge sous la sustente 
AB contrepesée en A. 


enveloppe (®) ; nous appellerons cette chainette-support, chainette minimale et 
nous réserverons 4 son point caractéristique la désignation g sans indice. 


Nous rappelons que la longueur Aq est la plus grande longueur que peut 
avoir une corde de chainette, tracée suivant une direction donnée et dont 


Vordonnée canonique a une extrémité a une longueur imposée, soit A a; cette 


longueur Aq se calcule en fonction de l’ordonnée imposée Aa au moyen de 
la relation: Aa/Ag = t, owt est la valeur minima que peut atteindre dans une 
chainette le rapport entre l’ordonnée canonique de l’extrémité d’une corde et 
la longueur de celle-ci tracée parallélement A une direction fixe. 

La longueur d’une sustente A B, contrepesée de la quantité A a al’extrémité 
A est donc au plus égale a Aalt. 

Nous appellerons chainettes-supports tendues, avec la désignation (9’), 
toutes celles dont le point caractéristique p est compris entre y et qj, et 
chainettes-supports laches, avec la désignation (y"), toutes celles dont le point 
caractéristique p” est compris entre gj et @. 
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La chainette-support de la suspension minimale (y) coupe la suspension lége 

_ lache en un point K et la verticale de l’appui B en un point b”; cette méme 

verticale coupe la courbe enveloppe (®) en un point b’. 

. Sur chacune des chainettes-supports ainsi sélectionnées n’importe quel 
point C, compris entre A et la verticale du point B, peut étre un point de 

charge, dont la charge est mesurée par 1: la longueur de l’arc compris entre le 


point C et l’extrémité C’ d’une corde C’B’ équipollente a CB. 

En conclusion de toutes ces remarques, on se rend compte que les seuls 
points de l’espace ott une charge fixée 4 un cable porteur-tracteur peut se 
trouver en suspension sont ceux compris a l’intérieur de la zone délimitée par 


la caténe age supérieure AB, par l’arc Agi de la caténe lége inférieure AB, 


par Vare Yo ob’ de la courbe enveloppe (®) et par le vecteur vertical b’B. A 
Vintérieur de cette zone des points d’accrochage on distingue la zone en crois- 
sant, comprise entre les deux caténes léges, dont chaque point C ne dispose que 
d'une seule possibilité de suspension et partant que d’une seule charge. A cette 
zone est accolée celle figurée par le triangle curviligne K b” B et le quadrilatére 
— curviligne gj b’b” K a Vintérieur desquels on dispose pour chaque point de 
deux modes de suspension et de deux charges. Les deux modules des caténes 
d’un point a double possibilité de charge sont inférieurs a celui, m, de la 
suspension minimale, si le point est dans le quadrilatére curviligne, tandis que 
dans le triangle curviligne Kb" B lun des modules est supérieur a celui, m, de 
cette suspension minimale (q), l'autre lui étant inférieur. 


2. Détermination de la valeur d’une charge en fonction du point de 
l’espace ou se trouve son point d’accrochage 


Pour explorer tous les points de la zone ot des charges peuvent étre accro- 
chées, il suffira d’examiner la valeur des charges pouvant étre fixées a l’extré- 


mité C d'une caténe de suspension AC tracée sur une chainette-support quel- 
conque de l’éventail délimité par les deux chainettes isosthénes (g) et (yo), 


issues de A et passant par B. 
Tout point C d’une telle chainette, compris entre A et la verticale de 


lappui B, détermine un vecteur rectiligne CB, auquel correspond dans la 

chainette considérée une, et une seule, corde équipollente C’B’. = 
La charge suspendable en ce point C est évidemment mesurée par l’arc CC’, 

dont l’extrémité C’ est celle de la corde C’B’, équipollente au vecteur CB. 


x 


Réciproquement a tout point B’ de la chainette, donc a un vecteur de 
translation B’B, ne peut correspondre qu’un seul point de charge C, marquant 
lextrémité de la caténe de suspension AC supportée par ladite chainette. 
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En effet si ce point B’ avait pu engendrer deux points de de charge Ce et Ce 
c’est qu’on aurait pu tracer 4 partir de B’ deux cordes B’ B'Ci et Bik respec- 
tivement équipollentes aux vecteurs BC, et poe ‘Mmarquant sur la chainette 
une corde, C,C iC», €quipollente 4 une autre corde, ee C3, ce qui est impossible. 

Quand le point générateur B’ se déplace de la droite (alors a l’infini) vers la 
gauche, son point de charge C se déplace également vers la gauche sur la 
chainette-support. S’il n’en était pas ainsi c’est que C pourrait rétrograder, donc 
étre engendré par deux points différents Bo et B', auxquels auraient corres- 


pondu sur la chainette deux cordes, BC et B,C, chacune équipollente au 


vecteur BC; ce qui est impossible. 
Mais il existe une position limite du point générateur B’, celle pour laquelle 
le point de charge C qu'il engendre, atteint I’ appui 4 A. Ce point limite est donc 


Vextrémité By de la corde Aj B; équipollente a AB. Dans ce cas la caténe de 


suspension AC se réduit A l’élément rectiligne infiniment petit, ayant pour 
direction celle de la tangente en A a la chainette-support étudiée, et la charge 


suspendue est égale au poids d’une longueur de cable égale a celle de l’arc AA}. 
On remarquera que le point limite B;, générateur du point de charge 
confondu avec le point A, ne peut étre qu’a droite de la verticale du point B. 
Dans le cas exceptionnel d’un point capable de deux suspensions, donc de 
deux charges (points situés a l’extérieur de la caténe lége lache, dans le triangle | 
semi-curviligne qo b’ B) il est facile de démontrer que la plus grande des deux 
charges correspond a la suspension ayant le plus grand module. 
En effet on devra tracer sur chacune des deux chainettes isosthénes (q,) et 


(~2) se coupant en C une corde équipollente au vecteur CB; soit Ci By celle 


tracée sur (p,), et CBs, celle tracée sur (q2). Si le module m, est plus petit que 
mz, Vordonnée canonique du point C; sera moindre que celle de C; cette 


inégalité, de méme sens que celle des modules, concourt a rendre l’arc CC , de 
(~,) plus petit que l’are CC} de (qq). 

Le cas particulier des charges suspendables en des points de l’espace infi- 
niment voisins des appuis se traitera facilement en utilisant les diagrammes des 


tractions exercées par un cable sur deux points fixes, A et B, auxquels il est 
attaché, quand on fait varier la longueur de ce cable’). 


nig M. NicoLas, Comment agit le cable tracteur d’un téléphévique, Z. angew. Math. Phys. 8, 
96-121 (1957). 
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3. Comment déterminer si une charge donnée 
peut étre accrochée au cable porteur-tracteur et alors 
en quels points de l’espace elle se trouvera 


Pour qu’une charge soit accrochable, il faut qu’il existe au moins une 
chainette-support dont un point C soit capable de cette charge. Or ce point C 
procéde d’un point générateur B’ sur la chainette, lequel détermine avec l’ex- 


trémité B de la sustente un vecteur BB’, auquel est équipollent une corde CC’, 


dont la caténe CC’ mesure la charge accrochable en C (voir figures 1 et 3). 


Figure 3 


Cas d’une chainette support, pour laquelle l’extrémité B est dans la zone monopode, Détermination 
du point a charge minima. 


Donc pour savoir si une charge, mesurée par le poids d’une certaine lon- 
gueur / de cable, peut avoir un point d’accrochage sur une chainette-support 
donnée (y’) ou (p”), de module m’ ou m", il suffira de tracer, dans le module de 
cette, chainette, une hypercirconférence de centre 6 et ayant pour hyperrayon 
cette longueur /, puis on examinera si elle a des points d’intersection avec la 
chainette. Ainsi grace 4 l’hypergéométrie, le probléme a traiter devient aussi 
simple que celui de l’intersection d’une droite (chainette ou hyperdroite) et 
d’une circonférence (hypercirconférence)*). 

Inversement, étant donné une chainette-support, il est possible de déter- 
miner les charges extrémes, soit maxima, soit minima, susceptibles d’étre 
accrochées en des points appropriés de celle-ci. 

Toutefois entre le probléme euclidien de l'intersection d’une circonférence 
avec une droite et celui de l’intersection d’une hypercirconférence avec une 


4) M. M. Nicoras, Géométrie de la chatnette ou hypergéoméirie, Z. angew. Math. Phys. 3; 
122-141 (1957). 
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hyperdroite, il existe une différence essentielle, qui constitue toute la difficulté © 
des problémes ot interviennent des cables. Alors que d’un point pris pour | 
centre on ne peut tracer qu’une seule circonférence tangente 4 une droite, on — 


peut d’un point pris pour centre tracer soit une seule hypercirconférence tan- 


gente a une hyperdroite (chainette) soit trois hypercirconférences, soit trés 


exceptionnellement deux, selon que le centre se trouve dans la zone monopode, 


ou dans la zone tripode de la chainette ou sur la ligne de démarcation entre les ~ 


deux zones. Et la complication de cette triple éventualité est aggravée par le 
fait que l’une des hypernormales présente une interprétation mécanique dif- 
férente de celle des deux autres. 


Nous désignerons par la notation (J") cette ligne de démarcation, qui n’est . 


autre que la courbe enveloppe des hypernormales a la chainette-support 
considérée, chaque point caractéristique y d’une hypernormale marquant le 
centre de l’hypercirconférence, qui au point ou elle est tangente a la chainette 
(point origine de l’hypernormale ou pode) y a méme rayon de courbure que 
celle-ci. 

L’étude du comportement d’un cable porteur-tracteur consiste donc a 
distinguer les divers cas de tangence qui peuvent se produire entre l’extrémité 
B, non contrepesée de la sustente, et les diverses chainettes-supports comprises 
dans l’éventail yj — gj, déterminé par cette sustente. 

On sait que dans le cas d’un point B tripode, les hyperrayons des trois 
hypercirconférences tangentes a la chainette sont, l’un court, l’autre moyen, 
le troisiéme long: nous désignerons leurs points respectifs de tangence par les 
lettres 6, (pode court), 8, (pode moyen), et 6; (pode long). Chacun de ces points, 
dits podes du point B, marque l’origine d’une hypernormale (lieu des centres 
de toutes les hypercirconférences tangentes 4 une chainette en un point de 
celle-ci) passant par B (voir figure 4). 

On se rappellera aussi que les podes f,, et 8, d'un point B sont toujours par 
rapport a l’axe de la chainette du c6té opposé a celui ot se trouve ce point, 
tandis que le pode f, est toujours du méme cété (et méme toujours a une plus 
grande distance de l’axe que ne l’est le point B) et cela quelle que soit la nature 
podique du point B. 

Si le point B est dans la zone monopode, son pode unique sera encore 
désigné par la lettre 6, parce que, comme nous le verrons par la suite, l’inter- 
prétation mécanique de ce pode unique ne différe pas de celle du pode court 
d’un point tripode. 

Si le point B est sur la ligne des points bipodes (enveloppe des hypernor- 
males a la chainette) le pode court £, subsiste évidemment, tandis que l’autre 
pode, qui résulte de la conjonction des podes f,,, et B, (conjonction qui se 
produit quand le point B tend a sortir de la zone tripode), sera désigné par la 
lettre B,,;. Le pode £,,, d’un point bipode, B, est toujours par rapport a l’axe 
vertical de la chainette, du cdté opposé a celui du point B et A une distance 
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x 


de l’axe supérieure a celle de B (sauf si B est sur l’axe, auquel cas il dispose 
_ de deux podes courts, symétriques, et d’un pode long marqué par le sommet 
de la chainette). 
Avant d’interpréter les divers cas de podisme qui peuvent se présenter entre 
“une extrémité de sustente et ses diverses chainettes-supports, il convient 
d'interpréter, au sens du probléme étudié, les conséquences d’une tangence en 
un pode court, en un pode moyen, en un pode long. 


Figure 4 


Cas d’un point B en zone tripode d’une chainette support (p"). Ses 3 hypernormales de podes f,, B; et 
B,, de points caractéristiques y¢, V1 et Ym’ Ses 3 hyperparalléles de points de rebroussement Ro, R, et 
R,, (on n’a figuré que la moitié de chaque hyperparalléle). 


Le cas du pode court est le plus simple: tout allongement de Pay pemayon 
fait apparaitre sur la chainette et de part et d’autre de f, un point dinter- 
section avec la nouvelle hypercirconférence. Au contraire tout raccourcissement 
de l’hyperrayon de contact en f, fait disparaitre Vintersection, Donc Vhyper- 
rayon de contact au pode f, mesure la plus petite charge susceptible d’étre 
accrochée en un point de la caténe de suspension, sous réserve que ce point f, 
soit susceptible d’étre effectivement un point générateur de charge (nous avons 


ZAMP IXa/22 
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vu que pour qu'il en soit ainsi il faut que f, soit a droite de l’extrémité B;, de 
la corde A/B; équipollente a AB). 

En particulier le pode court d’une chainette-support, dont le sommet est 
4 droite de l’aplomb de l’extrémité B de la sustente, ne peut étre générateur de 
charge puisque, étant a gauche de la verticale de B, il est nécessairement a 
gauche de l’extrémité B; de la corde A; B; équipollente a AB. 

Donc si le pode court peut étre générateur d’un point de charge C, (a 
charge minima alors), on en déduit que de part et d’autre de celui-ci les charges 
accrochables vont croissant. 

Notre réserve quant a la possibilité d’obtenir sur une chainette-support un 
point de charge effectif, 4 partir du pode court A, de l’extrémité B sur cette 
chainette, nous conduit 4 remarquer que quand #, est susceptible d’étre pris en 
considération comme point effectivement générateur, les deux autres podes f,, 
et 8, sont automatiquement hors de compétition et vice versa (puisque quand 


a 


B, est a droite de l’aplomb de B, f,, et 6, sont a gauche, et vice versa). Nous _ 


remarquons aussi que quand les podes f,, et 6, sont susceptibles d’étre généra- 
teurs de points de charge effectifs, l’efficacité de 6, implique celle de f,,, (puisque 
6, étant alors a droite de l’extrémité B/ de la corde A; B/, équipollente a AB, le 
point £,, l’est aussi, puisque sur la chainette il est encore plus a droite que f;). 


Les podes, moyen f,, et long f,, s'interprétent en étudiant la progression le — 


long de la chainette des points d’intersection de celle-ci avec les hypercircon- 
férences de centre B au fur et 4 mesure qu’on augmente l’hyperrayon qui a 
produit la tangence au pode. 

Au début de la croissance de cet hyperrayon on n’a que deux points d’inter- 
section, qui vont en s’éloignant de #,, l'un numéro 1, s’élevant sur la branche 
de la chainette, l'autre numéro 2, tendant vers le sommet. Quand la croissance 
de cet hyperrayon lui fait atteindre la valeur du contact en f,,, apparait alors 
entre la chainette et l’hypercirconférence un troisiéme point commun, celui du 
contact en f,,, tandis que des deux autres points d’intersection du début, l’un 
(1) se sera encore plus éloigné de f,, tandis que l’autre (2) se sera rapproché 
de f;. 

Continuant a faire croitre l’hyperrayon, le troisisme point commun apparu au 
contact en f,, se dédouble en deux nouveaux points d’intersection, numéros 3 
et 4, générateurs éventuels de deux points de charge a charges croissantes de 
part et d’autre du point de charge C,,,, engendré par £,,. Continuant a faire croitre 
les hyperrayons, les points d’intersection, numéros 1 et 4, s’éloignent respec- 
tivement vers le haut sur leur branche de chainette, tandis que les points 
dintersection 2 et 3 se rapprochent et tendent vers le pode f;. Donc celui-ci 
est générateur éventuel d’un point de charge C,, capable de la plus grande de 
toutes les charges susceptibles d’étre accrochées en un point de la chainette- 
support considérée; de part et d’autre de ce point C, les charges accrochables 
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_ sont évidemment égales deux a deux, tout en diminuant au fur et 4 mesure 
qu'on s’éloigne de C,. 

Ainsi quelle que soit la position relative de l’extrémité B de la sustente a 
_ légard d’une chainette-support de 1’éventail (@%) — (y/’), on est toujours str 
_d’avoir a considérer sur celle-ci un point générateur éventuel d’un point a charge 
minima, de part et d’autre duquel les charges accrochables égales deux a deux 
vont en augmentant: ce point remarquable, est, selon les cas, engendré soit 
par le pode f, (et alors la forme de la chainette-support est telle que son sommet 
est a gauche de l’aplomb de I’extrémité B), soit par le pode £,, (et alors le 
_ sommet de la chainette-support est a droite de l’aplomb de l’extrémité). 


x 


4. Ligne des points 4 charge minima 


Ainsi chaque chainette-support (p,) comporte un point éventuel a charge 
minima, C*,, de part et d’autre duquel les charges accrochables vont en aug- 
“mentant. I] est évident que ces points a charge minima forment dans le plan 
_ vertical de la sustente une ligne continue — que nous désignerons par la notation 
(QV ) —joignant le point C;? de la caténe lége (py) a celui C7? de la caténe lége (v7), 
et il est non moins évident qu’en ces deux extrémités les charges minima 
suspendables sont nulles, ce qui fait présumer que de C/? a C,’° les diverses 
charges minima afférentes aux chainettes-supports de l’éventail gp — gy vont 
d’abord en croissant pour atteindre sur une certaine chainette-support leur 
plus grande valeur, en un point Cy, puis décroissent jusqu’a la valeur nulle. 
C’est ce que nous démontrerons par la suite. 

Anticipant sur cette démonstration, tirons en immédiatement les consé- 
quences pratiques. Pour qu’a une charge donnée on puisse faire correspondre 
dans l’espace une ligne continue de A a B (cheminement statique) de ses points 
d’accrochage au cable supposé a l’arrét, il faut nécessairement que cette charge 
soit inférieure 4 la plus grande des charges minima, sinon le cheminement ne 
pourrait franchir la ligne de démarcation. 

On peut donc comparer a un filtre, a l’égard des charges, le réle de cette 
ligne de démarcation. 

La chainette-support (~j,), dont la charge minima accrochable (Cj,) est la 
plus grande de toutes celles minima accrochables aux diverses chainettes- 
supports, divise le champ des points de charge possibles en deux zones, supé- 
rieure et inférieure, ou l’on trouve toujours accouplées deux chainettes-sup- 
ports (y,) et (y,), l'une en zone supérieure, l’autre en zone inférieure, dont les 
points 4 charge minima respectifs c* et Ci, supportent des charges égales (et 
évidemment inférieures a (C,)). 

Nous verrons par la suite que cette chainette-support remarquable (p,,) 


est nécessairement une chainette-support lache. 
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A toute charge minima commune a deux chainettes-supports (px) et (5) 
correspondront deux cheminements statiques, qui, sauf aux points de fran- 
chissement C*, et Ci, de la ligne de démarcation (Q/), seront entiérement tracés 


lun au-dessus de la chainette-support (,) de la zone supérieure et l'autre au- — 
dessous de la chainette-support (;) de la zone inférieure. Donc le cheminement — 
statique inférieur d’une charge pouvant franchir la ligne (QW) des minima 
n’atteindra jamais l’appui B non contrepesé, et les seuls cheminements statiques — 


joignant A a B sont ceux franchissant la ligne (QC) au-dessus du point Cy, 
| 
| 
| 
| 
| 


Figure 5 
Forme présumée de la ligne (Ofl) — cas d’une sustente < AB,,. 


c’est-a-dire composé de points de charges situés sur des chainettes-supports 
d’un module supérieur a celui de la chainette (pf,). On peut méme préciser que 
le long de ce cheminement statique les chainettes-supports mises a contribution 
ne le sont que deux fois, pour des points de charge sis de part et d’autre de la 
ligne (OZ) (voir figure 5). 

Ce bref apergu des conséquences mécaniques de l’existence de cette ligne 
(Q) montre Vutilité d’une connaissance approfondie de celle-ci pour traiter 
tout probléme ot interviennent des cables porteurs-tracteurs. Aussi convient-il 
d’en exposer au moins une étude sommaire. 

Nous commencerons par mettre en place les points extrémes C/? et €//°, 
afférents aux chainettes-supports (qj) et (pi), des deux caténes léges AB. 


Pour chacune de ces deux chainettes-supports extrémes (yj) et (2), ou plus 
précisément pour des chainettes-supports de I’éventail g — qj, qui leur sont 


extrémité, Cc’? 
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infiniment voisines, l’extrémité B est dans la zone monopode (et méme dans 
la partie de cette zone qui est 4 droite de l’axe de la chainette). 

D’autre part, puisque le tracé de toute hypernormale, au voisinage de son 
origine, est normal a la chainette, on en déduit que le vecteur BBo, qui joint 
lextrémité B a son pode #, sur l’une ou l’autre des chainettes-supports 
extrémes (gy) et (yo), est un vecteur infiniment petit, dont la direction est 
normale a celle des tangentes en B, respectivement a la chainette (yj) et a la 
chainette (7). 

Nous remémorant alors la détermination du point de charge engendré par 
un point B’ (en Voccurrence f,, confondu avec B), nous en déduisons que 
Vextrémité C,/° de la ligne (QZ) est marquée par le point de la chainette (yj), ou 
la tangente est orthogonale a celle en B a cette méme chainette, et que l’autre 

m» est marquée par le point de la chainette (yj), ot la tangente 
est orthogonale a celle en B a cette méme chainette. 

Nous remarquons que le point C,”° a toujours une existence effective en 
tant que point de charge, puisque toujours situé a droite de l’appui A, sur la 


caténe lége lache AB. (En effet l’angle polaire — celui des deux tangentes a 
Yextrémité d’une caténe — d’une caténe lache étant nécessairement aigu, le 
point C,/° ne peut étre au-dela de A.) 

Au contraire le point C/? n’aura une existence effective, en tant que point 
de charge, que s’il est en deca de A. Or l’angle polaire d’une caténe tendue peut 
fort bien étre orthogonal et méme obtus. (L’étude des caténes minimales 
apprend que leur angle polaire est toujours aigu, sauf dans le cas extréme de la 
corde minimale verticale, et l’étude des caténes orthogonales apprend que 
celles-ci ont toujours leur ple au-dessus de la base de la chainette; donc parmi 
les caténes tendues dont la corde est paralléle a une direction donnée, figure 
toujours une, et une seule, qui soit orthogonale.) 

Ainsi selon que la caténe lége tendue AB aura un angle polaire aigu ou un 
angle polaire obtus, l’extrémité C/? de la ligne (Q//) des points a charge minima 
sera en de¢a (et alors C;? est un point de charge effectif) ou au dela de A. 

L’étude de la ligne (07) apprend aussi que quand son extrémité supérieure 
C/o est au dela de A, la ligne (7) passe alors par A; autrement dit parmi les 
chainettes-supports de l’éventail y — go, correspondant a ce cas de sustente, 
ce n’est qu’a partir, et en dessous d’un certain module, que le point de charge 
minimum d’une chainette-support se trouve a droite de A, et est alors efficace 
en tant que générateur d’un point de charge utilisable. ce 

Si nous reprenons la représentation de la sustente AB, contrepesée de A a, 
et dont la plus grande longueur est Ay = Aa/t, on voit que parmi toutes les 
extrémités -possibles entre A et y il s’en trouvera une, B,, qui donnera une 


caténe lége tendue AB,, orthogonale. 
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Cette longueur AB,, ou plus exactement le rapport AB,,| Aa= t,, se déter- : 
mine sans aucune difficulté en fonction de l’angle que fait la sustente avec © 


Vhorizontale. Selon que cet angle est positif (sustente ascendante) ou négatif 


(sustente plongeante), 4 valeur absolue égale de cet angle, la longueur AB,,, 


ascendante, est plus grande que la valeur AB,, plongeante. (Par exemple pour ~ 


une sustente horizontale t,, = 1,246 et pour une sustente inclinée 4 15°15’30" le 
méme rapport vaut 1,520 en sustente ascendante et 1,036 en sustente plon- 
geante.) 

Dans le cas d’une sustente dont la longueur est moindre que la longueur 


AB,, susvisée (c’est-a-dire dont la caténe lége tendue AB a un angle polaire 
aigu), il est facile de déterminer la chainette-support dont le point a charge 
minima coincide avec l’appui A, et en dessous de laquelle (donc pour des 
modules moindres) les points a charge minima se trouvent en deca de l’appui A. 
Sur cette chainette-support l’extrémité B; de la corde A; B;, équipollente a la 


sustente AB, devra étre le pode du point B, c’est-a-dire l’origine d’une hyper- 
normale passant par B. 
Or la construction d’un point quelconque d’une hypernormale dont on 


connait l’origine J sur une chainette se fait en associant deux points, A et H, de 
la chainette qui soient, par les longueurs des caténes, conjugués harmoniques des 


extrémités de la caténe orthogonale I] (Z point de la chainette dont la tangente — 


est orthogonale a la direction de la tangente en J); deux points, A et H, 
conjugués harmoniques d’un couple orthogonal J (un seul point suffit pour le 
définir) satisfont donc a la relation Al/AJ = — H1/ HJ. De ces deux points 
lun, A, est nécessairement au-dessus de la corde orthogonale LJ, l’autre en 
dessous mais sans cependant dépasser le milieu K de la caténe ae partir de 
A, on trace alors un vecteur AB, équipollent a H/, et l’on obtient un point B 
de I’hypernormale du point J; ainsi se trouve réalisée A Végard de cette 
chainette une sustente AB pour laquelle la corde équipollente aboutit bien au 
pode du point B (voir figure 6). 

On voit qu’a tout point J’ d’une chainette quelconque, prise comme chai- 
nette de référence, on peut faire correspondre une corde, dt ayant une 
direction imposée, en l’occurrence paralléle a celle de la sustente, sous la seule 
réserve que le point / ¢ soit tel que le milieu A’ de la caténe orthogonale RP 


détermine une corde K’ J’ faisant avec V’horizontale un angle supérieur a 
celui de la sustente. La détermination de ce point J’, limite, se fait sans 
difficulté, en utilisant notamment cette particularité géométrique des caténes 
orthogonales dont la longueur est toujours égale au double de l’ordonnée 


canonique de leurs milieux K’ soit: I’ J’ = 2 K’k’. 
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Ceci étant, on se rend compte facilement de ce que pour chaque point J’, 


pode présumé d’une extrémité B’, la corde H’ J’, paralléle a une direction fixe, 
sera d’autant plus petite que J’ sera plus bas sur la chainette, tandis qu’au 
contraire le conjugué harmonique 4A’ de ce point H’ s’éléve; donc plus le point 
J est bas, plus le rapport H’J’/A’ a’=A’ B’/A’ a’ diminue. 

Pour transposer dans le faisceau des chainettes isosthénes la combinaison 
réalisée sur cette chainette de référence, il faut en modifier le module dans le 


Figure 6 


Construction point a point d’une hypernormale a partir de conjugués harmoniques A, H du couple 
orthogonal IJ. 


rapport AalA’ a’. Donc plus sera courte la sustente AB (contrepesée de A a) 
et plus sera lache (module plus petit) la chainette-support pour laquelle le pode 
de l’extrémité B engendre un point de charge confondu avec l’appui A. Et 
pour cette sustente, seules les chainettes-supports d’un module moindre que 
celui susdéfini auront des points 4 charge minima en deca de l’appui A. 

En particulier quand le point J’ détermine une caténe orthogonale 1’ /’ 


dont le milieu K’ détermine une corde, Ki gh: paralléle a la direction fixe (celle 
de la sustente), le point H’ étant confondu avec K’, son conjugué harmonique A 


. . . A f b / y "4 
est éloigné a l’infini sur la chainette et son ordonnée canonique A’a’ étant 
devenue infiniment grande, le rapport H’J’/A’a’ tend vers zéro, donc la 


sustente AB devient infiniment petite. 
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Ainsi se trouve indirectement démontré le fait que la ligne de démarcation 


(QC) passe par l’appui A, quand la caténe lége tendue AB aun angle polaire 
obtus. 


5. Détermination de la chainette-support dont la charge minima 
est la plus grande de toutes les autres charges minima 


Il nous reste maintenant a étudier la détermination de la chainette-support, 
dont le point de charge minima supporte la plus grande de toutes les charges 
minima afférentes aux diverses chainettes-supports de l’éventail. 

Cette chainette-support est celle dont le point caractéristique gy” marque 
sur la courbe enveloppe (®) des chainettes-supports le point dont l’hypernor- 
male isomodulée a cette courbe (®) (c’est-a-dire tracée suivant le module m” 
de la chainette dont le point caractéristique est py”) passe par l’extrémité B 
de la sustente. 

La démonstration de cette propriété fondamentale des cables porteurs- 
tracteurs se déduit trés simplement de l'étude des hyperparalléles a la courbe 
enveloppe (). 

Mais faute de pouvoir faire dans le cadre de cette étude celle compléte des 
hyperparalléles a la courbe enveloppe (®), nous proposerons une autre démons- 
tration qui ne recourt pas, ou que trés peu, a l’hypergéométrie. 

Cette démonstration est fondée sur la connaissance des propriétés d’une 
certaine courbe, dérivée de la chainette suivant le processus décrit ci-aprés. A 
tout point A d’une chainette, on fait correspondre un point, B, extrémité d’un 


vecteur AB paralléle a une direction fixe, et dont la longueur AB est dans un 


rapport constant, AB/Aa, avec l’ordonnée canonique Aa du point A. Nous 
appellerons cette courbe la transformée oblique de la chainette, suivant telle 
direction, et dans tel rapport, chacun de ses points B, étant le transposé 
oblique d’un point, A, de la chainette, moyennant un vecteur de transposition, 
AB, paralléle une direction fixe et orientée et ayant une longueur proportion- 
nelle 4 l’ordonnée canonique de A (voir figure 7). 

L’étude des transformées obliques est trés féconde et si facile, qu'il nous 
parait suffisant d’en citer leur propriété fondamentale: en tout point B d’une 
transformée oblique, la tangente A cette courbe coupe la base la ou celle-ci est 
coupée par la tangente a la chainette au point A, dont B est le transposé 
oblique. 


Soit donc AB une sustente contrepesée de la quantité Aa, soit Ap la 
longueur maxima que peut atteindre une telle sustente ; prenons comme 
chainette de référence la chainette pour laquelle la sustente Ag, associée a 
Vordonnée canonique A a, marque aussi la corde minimale (les tangentes aux 


extrémités de la caténe Ap se coupent en T sur la base) 
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_ Tragons a partir de cette chainette de référence, de module m, la trans- 
formée oblique (8) “suivant la direction de la sustente (donc celle de Aq) et 
dans le rapport A B/A a, qui définit la longueur AB de la sustente en fonction 


de sa contrepesée A a. 
Cette courbe (.“4) coupe nécessairement la chainette en deux points, puisque 


AB/Aa < Ag/Aa. En ces deux points d’intersection, soit By au-dessus de 


Bi 


Figure 7 
Chainette et transformée oblique de chainette. 


Ag et Bi au-dessous de Aq, les cordes Al Bi et Aj, BY sont appareillées (c’est-a- 
dire que Aj By/Ag ay = Ao Bo/ Ap aj); de méme tout vecteur de transposition 
A"B" tracé au-dessus de ae mais, prolongé jusqu’a la chainette en un point 
By devient une corde A" B"; a tout vecteur A’ B’ situé en dessous de Ap on 
peut toujours associer un vecteur, et un seul, A” B” situé au-dessus de Ag, et 
tel que A’ B,/A’a’ = A" By/A" a"; il suffit en effet que les deux cordes A’ Bi, et 


A’ An By soient appareillées. 
Bats une telle transformation oblique un vecteur quelconque de transpo- 


sition A" B" (ou A’ {’ B’) associé avec la chainette de référence représente, a 
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Vhomothétie prés, une chainette-support de l’éventail gy — go, dont le module | 
est soit m= mx Aal/A"a" soit m’ =m x Aa/A’a’; les chainettes-supports | 


x N aries A " 
extrémes (ou les deux caténes léges) sont les homothétiques des caténes Ag Bo | 


et A, By de la chainette de référence. | 

De méme deux chainettes-supports (p’) et (p”) ayant un point commun — 
entre B et w sur le vecteur Ap sont les homothétiques des deux caténes décou-_ 
pées par deux cordes appareillées sur la chainette.de référence. 


Ainsi 4 chaque chainette-support de l’éventail @j — go d'une sustente AB, 


contrepesée de A a, correspond un, et un seul, vecteur 4” B” (ou A’ B’) qui est _ 
son vecteur figuratif sur la chainette de référence. 

A partir d’un point B" de la transformée (4) tracons l’hypernormale a la 
chainette de référence, qui aboutit 4 celui des podes 6” de B” qui est le plus a 
droite : nous savons que, quelle que soit la position du point B” (ou B’) par rapport 
aux zones monopode et tripode de la chainette de référence, ce pode B” a l’ex- 
tréme droite est soit un pode court du type f,, soit un pode moyen du type f,,. 


Pour ce pode l’hyperrayon de l’hypercirconférence de centre B” est la 
longueur J” de la caténe découpée sur la chainette de référence par une corde 


équipollente au vecteur B” 8”. Supposons qu’au vecteur 4" B" on en substitue 


un autre, infiniment voisin, A/ B; plus proche de Ag, donc ayant un module | 


plus grand: les deux points Aj’ et By! se trouveront sur les tangentes A” T” et 


B"T", tandis que le pode B/ du point B/! se trouvera sur la tangente B’@", ala 
tangente a la chainette au point 6”. Le nouvel hyperrayon J,’ correspondra donc 


a une corde équipollente au vecteur Bj B;. Si le point O” est a droite du point 


T", au vecteur B,’ 6,’ sensiblement paralléle au vecteur B” 8” correspondra une 
nouvelle corde équipollente qui pourra étre soit plus grande, soit plus petite que 
la précédente mais alors dans ce cas avec un rapport de réduction moins élevé 
que n’est le rapport d’accroissement des modules (soit m///m" = A" a"/A/ aj’) 
des chainettes dont les vecteurs figuratifs sont A" B" et Aj! By’. 

Donc, quand par les homothéties Aa/A” a" et Aa/A/’ aj’, on aura réalisé les 
deux chainettes-supports de modules respectivement égaux A mx AalA" a" et 
a mx Aa/Aj ax, Vhyperrayon minimum de la chainette-support m”k sera 
supérieur a celui de m", parce que quand on fait varier dans les mémes propor- 
tions la corde d’une caténe et le module de celle-ci, le rapport caténe/corde de 
leurs longueurs demeure inchangé, tandis que si la corde s’allonge plus vite 


(ou diminue moins vite) que ne croft (ou décroit) le module de la caténe, le 
rapport caténe/corde augmente. 


Ainsi de deux chainettes-supports infiniment voisines, figurées par les 
vecteurs infiniment voisins A” B" et Aj Bj, c’est celle de plus grand module 
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_ dont la charge minima est la plus grande, @ condition que la tangente a la chai- 
nette-support, de module m", au pode B" de l’extrémité B" de la sustente coupe la 
base en un point O" a droite du point T" ow celle-ci est coupée par la tangente a la 
chainette au point A". 

Un raisonnement identique, appliqué au cas ot le point @” est a gauche de 
T”, montre évidemment que les conséquences sont aussi inverses. On en déduit 

_ que la plus grande des charges minima est celle afférente a la chainette-support 

_ sur laquelle le pode f de l’extrémité B de la sustente marque le point ot la 
tangente coupe la base 1a ot celle-ci est coupée par la tangente en A a cette 

_méme chainette. Donc f est le conjugué minimal du point A ; donc la chainette- 

support, capable de la plus grande des charges minima, est celle dont le point 
caractéristique gy” marque le pode de l’hypernormale dismodulée a la courbe 

_enveloppe (®), qui passe par l’extrémité B de la sustente. On est stir que cette 
chainette-support fait partie de celles que nous avons qualifiées de laches, 

_ puisque le point 6”, pode du point B”, est nécessairement en dessous de la 


droite A” B” quand le point A” devient le conjugué minimal de ce point £”. 

On est sir aussi que cette chainette-support exceptionnelle est unique (et 
n'est nullement subordonnée a un choix du pode f a l’extréme droite), car, par 
un raisonnement fondé exclusivement sur l’hypergéométrie (en l’occurrence, 
l’étude de la courbe enveloppe des hypernormales dismodulées), on montre que 
par un point B quelconque, situé a droite de la verticale du nceud A d’un 
faisceau de chainettes isosthénes, il ne peut passer qu’une seule hypernormale 
dismodulée ayant pour pode un point situé sur la branche droite de la courbe 
enveloppe (®). [Une hypernormale dismodulée a la courbe enveloppe (®), est 
celle tracée a partir d’un point @ de celle-ci sous le module m de la chainette 

_isosthéne dont ce point @ est le point caractéristique. | 

D’ailleurs cette longue démonstration se réduit 4 un raisonnement exprimé 
en deux ou trois phrases quand on résout le probléme par l’hypergéométrie. 

Désormais nous distinguerons ce point unique de la courbe (®), dont l’hyper- 
normale dismodulée passe par l’extrémité B de la sustente, par l’appellation 
«pode dismodulé de ]’extrémité B», tout en nous rappelant que cette hypernor- 
male est commune a la courbe enveloppe (®) et a la chainette isosthéne, dont 
ce pode marque le point caractéristique. 


6. A quelles conditions le pode dismodulé de l’extrémité B 
d’une sustente peut-il étre générateur effectif d’un point de charge? 


Cette question mérite d’étre résolue, car suivant que sur la chainette-sup-_ 
port (p”), définie par le pode (p") de l’extrémité B sur la courbe (®), l’extrémité 
B; de la corde A; B;, équipollente a AB, sera a gauche ou a droite du point yp”, 
le point de charge engendré par celui-ci sera ou ne sera pas un point de charge 
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effectif; et s’il n’est pas un point de charge effectif, tous les points de charge 
minima répartis sur la ligne de démarcation (DL ) entre A et C,//° sont des points 
dont la charge minima décroit quand on passe ‘d'une chainette-support 4 une 
autre de moindre module. Donc dans ce cas, a toute charge suspendable au 
départ de l’appui A correspondra sitrement un cheminement statique, de A a B, 
pourvu qu’on s’engage sur une chainette-support de module supérieur, sans 
point a charge minima. 

Nous allons démontrer qu’a chaque inclinaison de la sustente correspond 
une longueur de celle-ci (donc un certain rapport AB/Aa) en dessous de laquelle 
4 toute charge suspendable aux abords de l’appui A correspondra une ligne 
continue de points d’accrochage 4 l’arrét (cheminement statique) entre A et B. 

Rappelons nous d’abord que le point de charge engendré par un point quel- 
conque B’ d’une chainette-support [en l’occurrence le pode pm” de la chainette- 
support (p”)] est un point C, qui marque sur la chainette considérée l’extrémité 
C d'une corde CC’ équipollente au vecteur BB’ (ou, ici, Bop’ ie 

Rappelons nous également que par un point, B, il ne passe qu’une seule 
hypernormale dismodulée a la courbe enveloppe (®). 

Enfin rappelons-nous également la construction, point a point, d’une hyper- 
normale, a partir des conjugués harmoniques A et H, par rapport au couple de 
points conjugués orthogonaux, dont l’un est le pode de l’hypernormale (voir 
figures 6 et 8). 

Si dans chacune des constructions afférentes aux diverses hypernormales, 
par exemple dans la construction de l’hypernormale de pode J,, nous nous 
astreignons a n’utiliser qu’un point qui soit le conjugué minimal de J,, soit le 
point A, (et son conjugué harmonique H,), nous construisons un point Y, de 
Vhypernormale, lequel marque alors l’extrémité d’un vecteur AY équipollent 
a la corde H,J,. Ce point Y; est donc le point de l’hypernormale issu en quel- 
que sorte du conjugué minimal du pode de l’hypernormale. L’étude du vecteur 
H,J,; apprend que le pivotement de la corde minimale A, J, dans un certain 
sens, autour de l’extrémité A, _provoque le pivotement dans le méme sens de 
la corde conjugée harmonique H, J, et que celle-ci passe par toutes les directions, 
depuis celle de la verticale descendante HJ, quand J. J est au sommet de la chai- 
nette, jusqu’a celle de la verticale ascendante HJ, quand J est a Vinfini a 
droite sur la chainette. On notera aussi que la direction du vecteur HF , orienté 
de H vers J, est toujours plus voisine de la verticale ascendante que ne lest 
celle de la corde minimale A /, orientée également de A vers /. 

Il existe donc une direction de corde minimale, et une seule, pour laquelle 


le vecteur AY a une direction donnée, en l’occurrence celle de la sustente qui 


nous importe. Ce vecteur AW (ou plus précisément le rapport AWV|Aa= = ty) 
peut étre calculé suivant le procédé que nous indiquerons ci-aprés. 
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Ceci étant, marquons sur la courbe enveloppe (@) le point gy”, correspondant 


a la corde minimale Ag”, dont la conjuguée harmonique H”@” est paralléle a 
la sustente, de direction Ag. 


; Sur Aq, portons le vecteur AB” = H" Ww" nous obtenons ainsi une extré- 
mité de sustente, B", dont le pode dismodulé y" engendre, en un point con- 
fondu avec l’appui A, le point de charge C,,, ayant la plus grande des charges 


ee ee ee ee ae 
I; h 
Figure 8 


Construction du point Y;, d’une hypernormale de pode J;, 4 partir de conjugués harmoniques, dont 
lun, A;,, est le conjugué minimal de J,. 


minima accrochables a toutes les chainettes-supports de l’éventail go — gy, 
afférent a cette extrémité B”. En fait, comme nous le préciserons plus loin, ce 
point de charge exceptionnel peut selon les cas correspondre soit a la plus 
grande des charges minima accrochables, soit a la plus grande de toutes les 
charges accrochables. 

Considérons maintenant un autre point caractéristique, pj’, en dessous de 
@". Son hypernormale dismodulée coupe la ligne Ap en un point B;,, sis entre 
A et B" (deux hypernormales dismodulées issues de la branche droite de la 
courbe enveloppe (®) ne se coupent jamais au-dessus de cette branche) (voir 
figure 9). 

Marquons sur la chainette isosthéne (y/') la corde H/ p/, conjuguée harmo- 
nique de la corde minimale Agi’; cette corde Hip,’ a subi un pivotement dans 
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le méme sens que celui de la corde minimale Ag} par rapport a la ligne Ap", 


donc le vecteur A, YY, équipollent a H/g/, coupera l’hypernormale dismo- 

dulée (p/) en un point Y{ situé entre By et —y. ' f 
L’hyperrayon de l’hypercirconférence de centre By, tangente en g,, sous 

le module m/, a la fois 4 la courbe enveloppe (®) et a la chainette-support (p31), 


Pil ie a oh le eh 


A x a u a 
sera mesuré sur cette chainette par une caténe encadrant la caténe AHy (ceci 


Figure 9 


B" (ou By) extrémité de la sustente dont le point de charge engendré par son pode dismodulé gest 
confondu avec l’appui A. 


est la conséquence de la construction d’une hypernormale par points conjugués 


harmoniques), donc la corde équipollente 4 Bi py, tracée sur la chainette (p/), 
donnera bien un point de charge C,,, inutilisable, puisque situé au dela de 
Vappui A. 

Par un raisonnement analogue, on démontre qu’au contraire pour les 
extrémités By, sises entre B" et ~, le point de charge C,, est effectif. 

I] est donc extrémement important de connaitre la longueur d’une sustente 
contrepesée, en dessous de laquelle l’on sera assuré que pour toute charge 
suspendable a l’abord immédiat de l’appui contrepesé A on dispose d’un 
cheminement statique de A a B. 

Désormais nous appellerons longueur de démarcation pour une sustente la 
longueur définie par cette extrémité, que nous désignerons par la lettre By, 
dont le pode dismodulé pm" engendre un point 4 charge minima confondu avec 
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l'appui contrepesé A. Cette extrémité By est donc celle a partir de laquelle 
les podes dismodulés des autres extrémités sont effectivement genérateurs de 
points de charge, si ces extrémités sont a droite de By, et ne sont plus généra- 
teurs effectifs pour les extrémités 4 gauche de By. 


La détermination du rapport AB,/Aa = = t,, propre a chaque direction, se 
fait en dressant systématiquement des tables completes des caractéristiques 


des cordes minimales, notamment en calculant pour chacune d’elles, AJ, en 


grandeur et direction, la corde HJ, dont l’extrémité H est conjuguée harmo- 
nique de A, par rapport au couple orthogonal défini par l’extrémité /. 

L’établissement de ces tables ne présente aucune difficulté. 

En particulier dans le cas d’une sustente horizontale ce rapport ty a la 
valeur 0,29233; dans le cas d’une sustente ascendante faisant avec l’horizontale 
un angle de 15°34’ ce rapport est égal a 0,52, tandis qu’il devient égal a 0,14 
pour une sustente descendante, faisant sensiblement le méme angle (15°17’) 
en dessous de l’horizontale. 


7. Nature exacte du point de charge 
engendré par le pode dismodulé pm” d’une extrémité B de sustente 


Jusqu’a présent nous avons toujours considéré ce point de charge comme 

étant un point a charge minima; ceci n’est exact que pour autant que le pode 

_dismodulé du point B sur l’enveloppe (®) — pode unique en tant que pode situé 

sur la courbe (®) — ne correspond pas, pour ce méme point B, a un pode long, f;, 
sur la chainette-support définie par /. 

Considérons la figure qui nous a servi a démontrer, en utilisant la trans- 
formée oblique, que la plus grande des charges minima (plus exactement que 
le point de contact de l’hypercirconférence de centre B, ayant le plus grand 
hyperrayon), était sur la chainette-support dont le point caractéristique était 
le pode du point B. Si la transformée oblique ne traverse pas l’axe vertical de 
la chainette de référence — et il suffit pour cela que AB/Aa soit plus grand que 
0,55242/cos 7, en appelant 7 l’angle de la sustente avec l’horizontale — on sera 
certain que le pode dismodulé f correspond a un pode court du type f,. Les 
charges exceptionnelles afférentes aux diverses chainettes-supports ne peuvent 
alors étre que des charges minima. 

Ce n’est que dans le cas ot le rapport AB/Aa est suffisamment petit, que 
certaines chainettes-supports peuvent introduire l’extrémité B dans la partie 
gauche de leur zone tripode. 

Pour savoir si le pode dismodulé de l’extrémité B d’une sustente marque 
sur la chainette-support, dont le point caractéristique est ce pode dismodule, 
un pode moyen ou un pode long (auxquels cas on a affaire soit a un point 
générateur d’une charge minima, soit 4 un point générateur d'une charge 
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maxima), il faut connaitre la place sur l’hypernormale (8) du point caracté-_ 
ristique y de celle-ci. (Il est évident qu'il s’agit du point caractéristique de- 
V’hypernormale considérée comme faisant partie de la famille des hypernor- 
males isomodulées a la chainette-support (f) et non comme faisant partie de la 
famille des hypernormales dismodulées a la courbe enveloppe (®); ainsi une 
hypernormale dismodulée a, du fait de sa double appartenance a deux familles 
d’hypernormales, deux points caractéristiques y;, et yais, le point caractéristique 
dismodulé y,;, étant toujours plus éloigné du pode f que le point caractéristique 
isomodulé y;,.) 

La mise en place de y;, se déduit facilement de la propriété géométrique 
suivante des cordes minimales: 4 toute corde minimale AJ d’une chainette 
quelconque il ne correspond en Peck et en grandeur (notamment rapportée 


a ordonnée canonique A Aa) qu’un seul vecteur Ay joignant le conjugué minimal 
A de J au point caractéristique y de l’hypernormale a la chainette en J. Ce 


vecteur Ay, dont la direction pivote dans le méme sens que celui de la corde 
minimale, devient tour a tour paralléle a toutes les directions. Plus précisément, 


quand J se déplace de gauche a droite, A partir du sommet, le vecteur 4 y 
passe de la direction verticale descendante a la direction horizontale (atteinte 
pour l’abscisse canonique: 0,751 de J), puis a une direction verticale ascendante 
(atteinte pour l’abscisse canonique 1,0364 de J), puis le pivotement continuant, 
le point y passe a gauche de la verticale du point A : pour ces points y la direction 


des vecteurs Ay ne peut plus étre celle d’une sustente (nous rappelons que nous 


figurons a gauche l’extrémité contrepesée A de la sustente AB) (voir figure 10). 
Donc pour chaque corde minimale il sera possible de déterminer, en direction 


et en grandeur, ce vecteur Ay, et le rapport Ay/Aa= = t, achévera de déter- 


miner pour la direction en cause la longueur AB de la sustente contrepesée, 
au-dessus de laquelle le point a plus grande charge exceptionnelle est un pode 
moyen ou court, tandis que pour les sustentes d’une longueur moindre, le pode 
dismodulé de l’extrémité B marque sur la chainette-support correspondante 
un pode long. 


Ainsi pour une sustente de longueur supérieure a AB, (plus précisément, 


pour celles dont le rapport A B/Aa est supérieur 4 un certain taux, t,, propre 
a la direction orientée de la sustente) la charge exceptionnelle (celle du point de 
charge engendré par le pode dismodulé de l’extrémité) mesure la plus grande 
des charges minima accrochables aux diverses chainettes-supports susceptibles 
d’étre prises en considération. 

Pour une sustente de longueur moindre que AB, la charge exceptionnelle 
mesure la plus lourde des plus fortes charges accrochables aux diverses chai- 
nettes-supports susceptibles d’étre prises en considération. Dans ce dernier 
cas la plus grande des charges minima théoriquement accrochable n’a plus de 
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point d’accrochage effectif (il serait A gauche de l’appui A) et la plus grande 
_pratiquement accrochable des charges minima est alors celle de la chainette- 
support, dont le pode isomodulé [donc qui n’est plus situé sur la courbe (D) | 
engendre un point de charge confondu avec l’appui A. Cette chainette-support 
‘se détermine aisément en fonction de la chainette-support dismodulée de 


Vextrémité B,: soit AB la sustente en cause, le module de celle de ses chai- 
_nettes-supports qui engendre un point de charge confondu avec le point A sera 


Figure 10 
Courbe enveloppe (/") des hypernormales 4 une chainette (on n’a figuré que la branche gauche), 


alm y|m Point d’inflexion | alm y|m 
td | de (I): pea Tes 
A, — 1,385 2,122 i. Ag — 1,869 3,319 
IE 1,0364 1,586 Sia 0,27, dle 0,751 1,295 
Vy — 1,385 2,93 X 2,122] 5 Vo = 0593 3,217 
qT; — 0,25 0 sks = 2,764. | 15 — 0,82 0 


égal a celui, m,, de la chainette-support dismodulée de B,, multipliée par le 
rapport AB/AB,, soit: m= m, x AB/AB,,. 

Il est donc essentiel de connaitre pour chaque direction de sustente les 
emplacements relatifs des extrémités remarquables B, et B, (voir figure 11, 
les lieux de ces points). 

L’étude particuliére de ces deux extrémités apprend qu’il existe une direc- 
tion de sustente, et une seule, pour laquelle ces deux extrémités sont confon- 
dues: cette direction est celle d’une sustente ascendante faisant avec l’horizon- 
tale un angle dont la tangente a une valeur trés voisine de 13,5/100 avec un 
rapport t, = t, égal a 0,405 a cinq milliémes prés (en utilisant une table de 
fonctions hyperboliques a cing décimales, on ne peut encadrer cette direction 
qu’entre deux directions voisines successives, l'une de 12,77/100, l'autre de 


14,7/100). 
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Pour les sustentes orientées (a partir de l’extrémité contrepesée A) faisant 
avec l’horizontale un angle dont la tangente est supérieure a 13,5/100, l’extré- 
mité caractéristique B, est 4 gauche de l’extrémité B,,; tandis que pour des angles 
moindres (et notamment pour toutes les sustentes plongeantes) l’extrémité B, 
est a droite de l’extrémité B,,. 

Donc pour les sustentes ae faisant avec l’horizontale un angle infé- 


rieur a 13,5/100, et ayant une longueur (pour ces sustentes AB, “ AB.) 


(2) 


298d 


D 


135/100 


OMx4a 
Figure 11 
Faisceau de chainettes isosthénes: 1° Lieu des points caractéristiques isomodulés Visy des hyper- 


normales, dismodulées de podes qm, a la courbe enveloppe (®). 2° Lieu des points w, des hyper- 
normales dismodulées. 


moindre que AB, (ou AB/Aa < t,), la plus grande des charges bénéficiant 
d'une ligne continue de points d’accrochage entre les deux appuis A et B sera 
celle la plus grande pouvant étre accrochée a proximité de l’appui contrepesé A, 
tandis que pour des longueurs supérieures a AB, la plus grande des charges 
bénéficiant d’une ligne continue de points d’accrochage sera moindre que la 
plus grande de celles accrochables a proximité de l’appui A. 

Pour les sustentes orientées faisant avec l’horizontale un angle supérieur 


& 13,5/100 (pour ces sustentes A = AB.) et ayant une longueur moindre 


que AB,, la plus grande des charges bénéficiant d’une ligne continue de points 
d’accrochage sera encore celle la plus grande de celles accrochables 4 proximité 
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de l’appui contrepesé A, et cette particularité subsistera (différence avec le cas 


précédent), tant que la longueur AB ne dépassera pas la longueur AB,, 

_ partir de laquelle la plus grande des charges, a ligne continue de points d’ oy 
_chage, sera moindre que la plus lourde de celles accrochables a proximité de 
aa appui A. 


Conclusion 


Avant de tirer les conclusions de cette étude du cable porteur-tracteur, il 
_ convient de parfaire une observation exprimée au début et signalant une certaine 
ressemblance entre les problémes de cables porteurs-tracteurs a l’arrét et ceux 
des cables fixes d’autres installations (ponts suspendus par exemple). 
Si dans ces autres installations il peut se trouver des cables, sur lesquels 
-s’exerce une traction constante (contrepesée) en un de leurs points de contact 
sur un appui spécialement aménagé, sur lequel le cable peut glisser sans frot- 
_tement (théoriquement), cette derniére circonstance ne se retrouve plus sur 
_Vappui voisin (ou un des appuis voisins) oti le cable est fixe; ce cAble a une 
_extrémité fixe. Au contraire le cable porteur-tracteur est censé glisser sans 
-frottement sur Vappui non contrepesé, ot: 11 subit une traction appropriée a 
son état de charge, avec un changement de longueur, également approprié, 
notamment entre le point de charge et cet appui. 
Donc pour étudier l’accrochage d’une charge donnée en un point donné d’un 
cable contrepesé, mais fixé en un autre appui, Bb, on devra rechercher sur les 


x 


deux cheminements statiques afférents a cette charge (s’ils existent, quitte a 


modifier la contrepesée) les points C, dont la longueur de caténe CB, mesurée 
sous le module de la chainette-support du point C (c’est-a-dire, l’hyperdistance 


. CB, mesurée sous ce module particulier) est égale a la longueur résultant du 
point d’accrochage choisi sur le cable. 


Inversement, connaissant la longueur CB, on déterminera les charges 
pouvant étre effectivement accrochées en ce point en recherchant sur chaque 


-cheminement statique les points dont l’hyperdistance dismodulée CB (dis- 
modulée, puisque le module de leur chainette de mesure varie avec le point) 
est égale a la longueur imposée. 

On remarquera aussi que dans certains cas de sustente, il peut exister une 
charge, qui soit, non plus aux abords des appuis, la plus grande de toutes celles 
accrochables et pour laquelle on ne dispose que d'un seul point de charge 
(d’une seule longueur CB). 

Finalement cette étude du cable porteur-tracteur a l’arrét nous a fait 
reconnaitre l’existence sur chaque direction de sustente contrepesée de quatre 
longueurs remarquables de sustente, et de chacune nous avons donné Vinter- 


prétation mécanique. 
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Nous les rappelons briévement, avec leurs notations: > oy ; 
ie Ag, longueur maxima, définie en fonction de la contrepesée Aa par la 


relation Ag = Aax1 /x, au-dessus de laquelle il est impossible de sustendre un 
cable, et pour laquelle il n’existe qu’une seule forme de cable sustendu, exclu-_ 
sive de toute charge; ; 
2° AB,, (ou AB,,|/Aa = t,) longueur de la caténe lége tendue dont l’angle~ 
polaire est orthogonal et qui marque entre les sustentes la différence mécanique © 
suivante: pour chacune des sustentes de longueur moindre, les suspensions se 
divisent en deux catégories, savoir celles des chainettes a modules supérieurs, 
sur lesquelles les charges accrochables croissent constamment au fur et a- 
mesure que leurs points d’accrochage s’éloignent de l’appui contrepesé A, et 
celles des chainettes 4 modules inférieurs, sur chacune desquelles il existe un 
point, dont la charge est la plus petite de toutes celles pouvant y étre accrochées ; 
oe AB, (ou AB,|Aa = t,) longueur pour laquelle, selon que l’inclinaison 
de la sustente au-dessus de Vhorizontale est supérieure ou inférieure a 13,5/100 
environ, il y aura, accrochable entre les appuis, soit une charge maxima maxi- 
morum, soit une charge maxima minimorum; pour les sustentes ayant des 
longueurs différentes, la disposition de leurs charges accrochables exception-_ 
nelles varie selon que leurs longeurs dépassent ou ne dépassent pas la longueur 


caractéristique A B,,; 


4° AB, (ou AB,/Aa = t,), longueur dont la comparaison avec celle susvisée. 
AB,, (plus petite que A8, pour les sustentes ayant au-dessus de l’horizontale 


une inclinaison supérieure a 13,5/100, et plus grande que AB,, pour les sustentes 
dinclinaison moindre), modifie notablement la disposition des suspensions a 
charge maxima. 

La mise en évidence des conséquences mécaniques d’une modification de 
longueur d’une sustente (sa contrepesée demeurant inchangée) dans le cas d’un 
probléme de suspension, pourtant réduit a la plus grande simplicité, convainc 
que, quelle que soit la méthode employée pour calculer les cables, ce calcul ne 
saurait étre conduit commodément et efficacement si on ne dispose pas de 
tables donnant pour chaque direction, de degré en degré au moins, les quatre 
valeurs fondamentales évoquées ci-dessus. Ni l’établissement de telles tables, 
ni leur exploitation ne présentent de difficultés. 


Zusammenfassung 


Die Feststellung der mechanischen und geometrischen Merkmale eines Trag- und 
Zugseiles ist entweder auf empirischem Wege oder durch Berechnung zu erlangen. 
Im letzteren Falle kann den sich aus der Berechnung ergebenden Formeln keinerlei 
Vertrauen geschenkt werden, sofern diese Berechnung nicht ausgeht von einer ge- 


nauen Kenntnis der theoretischen Gegebenheiten, deren wechselseitige Beziehungen 
man in Form von Gleichungen ausgedriickt zu haben beansprucht. 


Vol. [Xa, 1958 Maschinelle Auswerterechnung bei der Flugabwehr 359 


Da der einfachste Fall derjenige einer einzigen Last ist, welche an einem Trag- und 
_ Zugseil im Zustand des Stillstandes aufgehangt ist, war es angebracht, eine griind- 
liche und vollstandige Analyse einer solchen Aufhangung vorzunehmen. Dies ist der 
Zweck dieses Aufsatzes. 
Es wird darin festgestellt, dass jeder Neigungswinkel, in welchem das Tragseil auf- 
-gespannt ist, gekennzeichnet ist durch 4 unterschiedliche Werte, welche jeweils das 
Verhaltnis zwischen der Lange des Tragseiles und seinem Gegengewicht bestimmen. 
_ Der Verfasser ist der Auffassung, die Berechnung von Tragseilen lasse sich nicht 
bequem und wirksam durchfiihren, wenn man nicht iiber Tabellen verfiigt, die 
_ ubrigens leicht aufgestellt werden kénnen, welche diese 4 Werte fiir die in Frage 
kommenden Neigungswinkel angeben. 
Ferner stellt man fest, dass die dargelegte Theorie sich auch auf feste Seile anderer 
Anlagen anwenden lasst, vor allem auf die Tragseile von Hangebriicken. 


(Regu: le 11 octobre 1957.) 


Maschinelle Auswerterechnung bei der Flugabwehr 


Von ERNST SCHULTZE, Bern?) 


1. Einleitung 


Die optische Erprobung [1]?) von Feuerleitgeraten fiir die terrestrische 
Fliegerabwehr verfolgt den Zweck, Angaben itiber die Genauigkeit solcher 
Gerate im Normalbetrieb zu gewinnen, ohne dass dabei geschossen zu werden 
braucht. Dabei wird einerseits das Ziel mit speziellen Theodoliten fortlaufend 
vermessen, und andererseits werden die an den Skalen des Rechengerates oder 
Kommandogerates erscheinenden Ergebnisse, namlich die ballistischen Koor- 
dinaten des voraussichtlichen Treffpunktes in regelmassigen Zeitabstanden 
fotografisch festgehalten und nicht an die Geschiitze weitergeleitet. 

Will man nur das Kommandogerat (Kgt.) allein priifen (unabhangig vom 
Radar), so lassen sich die kostspieligen Theodoliten entbehren. Dafiir muss 
man dann die Eingangswerte des Kgt. fotografieren. Die Kenntnis der Ein- 
gangs- und Ausgangswerte des Kgt. erméglicht allein schon die Beurteilung 
der dynamischen Rechengenauigkeit, vorausgesetzt, dass sich das Ziel gentigend 
gleichformig bewegte. Schliesslich lassen sich fiir diesen Zweck auch Zielflug- 
zeug und Radar entbehren, wenn man dafiir die Eingangswerte kiinstlich 
durch einen Simulator erzeugt. Die Auswertung solcher Versuche, das heisst 
die Berechnung der Fehler der Ausgangswerte bei gegebenen Eingangswerten, 
bietet, auch bei Beschrankung auf das Allernotwendigste, cine ganz betracht- 
liche Arbeit. Denn um ein zuverlassiges Bild von der Genauigkeit eines Ket. 
zu bekommen, miissen Fliige in verschiedenen Héhen, in verschiedenen Ent- 


1) Hasler AG, Werke fiir Tetephonie und Prazisionsmechanik. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 374. 
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fernungen und mit verschiedenen Geschwindigkeiten gepriift werden. Da bei 
jedem Flug viele Bilder gemacht werden, sind bald einmal einige tausend 


Treffpunkte zu berechnen. Es wurde deshalb der Versuch unternommen, eine 


solche Auswertung auf der ERMETH durchzufiihren (Elektronische Rechen- 
Maschine der ETH). Die vorliegende Arbeit soll iiber die dabei angewandten 
Rechenmethoden berichten. 


2. Grundlagen und Bezeichnungen 


2.1 Anordnung 


Man hat als Eingangswerte die Kugelkoordinaten des Zieles und als Aus- . 


gangswerte die ballistischen Koordinaten des zukiinftigen Treffpunktes. Will 
man das ganze Feuerleitgerat priifen, so verwendet man als Eingangswerte die 
aus der Theodolitmessung resultierenden Kugelkoordinaten; bei Priifung des 
Ket. allein verwendet man die ins Kgt. eingehenden Kugelkoordinaten. 
Eingangswerte und Ausgangswerte werden am einfachsten gleichzeitig foto- 
grafiert, und zwar in konstanten Zeitintervallen At. (Die Wahl des Zeitinter- 
valles unterliegt bei Anwesenheit eines Ladeverzuges k der Beschrankung, 
dass k/At eine ganze Zahl sein muss, siehe Abschnitt 4.4.) Die fotografierten 
Werte werden nach der Entwicklung der Filme abgelesen und auf Lochstreifen 
oder Lochkarten festgehalten. 
Die Auswerterechnung auf der ERMETH geht dann folgendermassen vorsich: 
1. Eingeben des auf Lochstreifen oder Lochkarten gespeicherten Rechen- 
programmes in die Rechenmaschine. 
2. Eingeben der auf Lochstreifen oder Lochkarten gespeicherten Eingangs- 
und Ausgangswerte eines Fluges in die Rechenmaschine. 
3. Durchfiihrung der Auswerterechnung, und zwar: 
a) Glattung der Eingangswerte, das heisst Berechnung des Normfluges 
(Abschnitt 2. 2), mit anschliessendem Drucken der Eingangsfehler; 
b) Schussfehlerberechnung auf Grund des Normfluges und der Ausgangs- 
werte (Abschnitt 2.3) mit anschliessendem Drucken der Schussfehler ; 
c) Drucken der Summenverteilungen der Eingangs- und Schussfehler. 
Dies gilt fiir den ersten Flug. Beim zweiten und den folgenden Fliigen 
braucht das Rechenprogramm nicht nochmals in die Maschine gegeben zu 
werden, sondern man kann dann gleich mit dem Eingeben der Ein- und Aus- 
gangswerte beginnen. 


2.2 Die Eingangswerte 


Da die Eingangswerte des Kgt. normalerweise mit Fehlern behaftet sind [2], 
miissen sie zuerst geglattet werden. Unter der hier gemachten Voraussetzung, 
dass das Kgt. nur fiir den gleichférmigen Flug gepriift werden soll (gleichf6r- 
miger Flug = Flug mit konstantem Geschwindigkeitsvektor), wird zundchst 
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auf Grund der Eingangswerte ein Normflug berechnet. Das ist derjenige gleich- 
férmige Flug, der gegeniiber den vorhandenen Eingangswerten die kleinsten 
Unterschiede aufweist. Er braucht bei Verwendung der Eingangswerte des 
Ket. nicht mit dem tatsachlich erfolgten Flug tibereinzustimmen (zum Beispiel 
‘wegen Radarfehlern), hingegen ist er fiir die Beurteilung der Rechengenauig- 
keit des Kgt. massgebend. Die Unterschiede zwischen den Eingangswerten des 
Kgt. und denen des Normfluges werden als Eingangsfehler bezeichnet. Der Ein- 
gangsfehler ist ein Vektor. Es kénnen von ihm sowohl der Absolutbetrag As, als 
auch die Komponenten berechnet werden. 


2.3 Die Ausgangswerte 


Wegen seiner eindeutigen Geschwindigkeit gestattet der Normflug eine 
eindeutige Treffpunktberechnung. Der Treffpunkt T fliegt dem Messpunkt M@ 
bekanntlich immer um die Summe von Geschossflugzeit t; und Ladeverzug k 
_voraus. Seine ballistischen Koordinaten «7, er und tr werden den Ausgangs- 
werten des Ket. gegeniibergestellt und die Differenzen als Kgt.-Fehler bezeich- 
net. Weiter bestimmen die Ausgangswerte a>, ¢; und t; des Kgt. eindeutig 
einen bestimmten Punkt im Raum, namlich den Sprengpunkt S (Ort der 
- Geschossexplosion), dessen ballistische Koordinaten sie sind. Der Abstand T'S 
der beiden Punkte ist der theoretische Schussfehler As, Als wirklichen Schuss- 
_fehler As, bezeichnet man hingegen den Abstand GS, wobei G der Gegnerpunkt 
ist (in welchem sich das Ziel zur Zeit der Geschossexplosion befindet [1]). Beide 
Schussfehler haben ihre praktische Bedeutung: Der theoretische Fehler ist 
identisch mit dem eigentlichen Kgt.-Fehler, der wirkliche Fehler ist die prak- 
tische Auswirkung des Kgt.-Fehlers. 

Es liesse sich noch ein dritter Schussfehler berechnen, der bei Verwendung 
von Annaherungsziindern als wirklicher Fehler in Frage kame. Dieser besteht 
im kiirzesten Abstand, den Ziel und Geschoss in irgendeinem Zeitmoment 
haben kénnen. 

2.4 Bezeichnungen 


Koordinaten: 
A Lagewinkel | 
e Distanz (Schragentfernung) Kugelkoordinaten, 
OL Seite 
er Kartenentfernung | Zylinderkoordinaten, 
4 mone : ‘ \ kartesische Koordinaten, 
x,  kartesische Grundrisskoordinaten J 
a Seite (inkl. Derivation) | 
€ Elevation | ballistische Koordinaten. 
G Geschossflugzeit 


362 


ERNST SCHULTZE ZAMP 


Alle diese Koordinaten (mit eventueller Ausnahme der ballistischen) be-— 
ziehen sich auf das Feuerleitgerat als Nullpunkt. Die x-Achse fallt mit der ~ 
Nullrichtung fiir die Seite zusammen (zum Beispiel nach Norden). Die Seite 
wird im Gegenuhrzeigersinn positiv gezahlt, das heisst, die y-Achse hat eine 


Seite von 1600 A°/oo. 


BOMHYZE 


Punkte: 


Messpunkt (Eingangswerte), _ 
Normpunkt, 

Treffpunkt, 

Sprengpunkt (Ausgangswerte), 
Gegnerpunkt, 

Wechselpunkt. 


Ein grosser Buchstabe, als Index an eine Koordinate oder einen Vektor 
angehangt, bezeichnet immer den Punkt, auf den sie oder er sich bezieht. Der 
Index 0 bedeutet, dass die betreffende Grésse zur Zeit ¢ = 0 gilt. 


7 
v 


Vektoren: 


Ortsvektor eines Punktes (x, y, h) , 
Geschwindigkeitsvektor mit den Komponenten (v,, v,, V,). 


Kenngrossen des gleichf6rmigen Fluges: 


Kurswinkel, 
Horizontalkomponente von v , 
Steigwinkel, 

Zeit des Passierens von W, 
Kartenentfernung von VW, 
Hohe von W. 


Fehler: 


Betrag des Eingangsfehlers, 

Komponente in Kursrichtung, 
Komponente senkrecht zur Kursrichtung, 
Komponente in Héhe, 


Fehler der Kugelkoordinaten, 


| Ket.-Fehler, 


theoretischer Schussfehler, 
wirklicher Schussfehier. 
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Sonstige Gréssen: 
t Zeit, 
At Bildintervall, 
k Ladeverzug, 
B = * — a = Flugrichtungswinkel (Winkel zwischen Visier 
und Kurs im Grundriss), 
Derivation, 


Einheitsvariable [Formel (4. 2)]. 


Ss m1 
——— 


3. Glattung und Eingangsfehler 
3.1 Berechnung des Normfluges 


Die Glattung, das heisst die Berechnung des Normfluges, wird am einfach- 
sten in kartesischen Koordinaten vorgenommen, denn in einem kartesischen 
System sind die Koordinaten eines gleichférmigen Fluges lineare Funktionen 
der Zeit: 

fy) = ty + at, Vy = vy Fat, hy Q=hy tut. 3.1) 


Die Berechnung der sechs Konstanten xy,, Vy,, Ay, Vx) Vy, Vy, die den Normflug 
vollstandig bestimmen, erfolgt durch Anwendung der Methode der kleinsten 
Fehlerquadrate ([3], § 21) gesondert auf jede der drei Koordinaten. 

Zur Umrechnung der Kugelkoordinaten auf kartesische dienen bekanntlich 


die folgenden Formeln: 


Crap Cx, COSAny sha = Sine COSKyy | (3.2) 


hy = Cy SMAyY, Vy = ey SINM, - | 


Die Kartenentfernung ¢,,, ist dabei eine Hilfsgrésse. Die Eingangswerte sind 
numeriert, wobei das erste Bild die Nummer 0 bekommt: 


Oy 5 


Ay; | Eingangswerte des Bildes Nr.7 (7 =0,1,2...,m”). (3.3) 


eu; 


Das Bild Nr. 7 ist zur Zeit 
t= 7 At (3.4) 


gemacht worden. Die kartesischen Koordinaten sind dann entsprechend: 


a a ae Ay, (PesOnd 2, vt. 5 7): (3.5) 


364 ERNST SCHULTZE ZAMP 


Nun bestimmt sich zum Beispiel die x-Komponente des Normfluges aus der 
Bedingung, dass 


n 


S = YF (xy (ts) — ty)? = Dy Cin, + Ye AtT — He,)? (3.6) 
j=0 


j7=0 
ein Minimum sein soll. Man hat also die beiden Gleichungen 


SOL EE Re eC (3.7) 


eee , = 
On Ni OV» 


welche durch Ausrechnen das folgende Gleichungssystem liefern: 


DP ky. Df ty AD ee De eT 0 At = DP i ae oe 


Daraus lassen sich xy, und v, bestimmen. Auf genau dieselbe Weise erhalt man 
auch die andern Komponenten y,y,, v, und hy,, v,, womit der Normflug voll- 
standig bestimmt ist. 

Die maschinelle Verwertung der abgeleiteten Formeln und iiberhaupt der 
ganzen Auswerterechnung erfolgt durch haufige Verwendung von Unterpro- 
grammen. Fiir jede gréssere Teilrechnung, die an mindestens zwei verschie- 
denen Stellen gebraucht wird, lohnt es sich, ein Unterprogramm (UP) herzu- 
stellen. Dadurch wird das Hauptprogramm wesentlich kiirzer und iibersicht- 
licher. Eine Zusammenstellung der verwendeten Programme befindet sich in 
Abschnitt 6. So erfordert die Berechnung des Normfluges die UP 18, 6 und 9. 
Dabei braucht UP 6 selbst zweimal das UP 4, und dieses wiederum zweimal 
das UP 2. 


3.2 Die Kenngréssen des Normfluges 


Die in den Gleichungen (3.1) auftretenden Konstanten Kn» VNy» Yn,» Veer Vy 
sind noch nicht die Kenngréssen des Normfluges. Als solche gelten vielmehr x, 
Up, Vy SOWIE ty, exw, hy. Die letzten drei fixieren die Lage des Fluges in bezug 
auf das Ket. (der Wechselpunkt W ist derjenige Punkt des Fluges, fiir den die 
horizontale Projektion den kleinsten Abstand vom Ket. hat), wahrend die 
ersten drei vom Kgt. unabhiangige Gréssen sind. 

Die Berechnung von x und v, ist einfach und kann wegen 


Vez = Uz COS%, Vy = VU, Sinx 


mit Hilfe von UP 5 geschehen. Die Berechnung von ty und egy geschieht an- 
hand der Figur 1 (Grundriss). Ny ist dabei derjenige Punkt des Normfluges, in 
dem sich das Ziel zur Zeit t = 0 befindet. 

Man hat die Formeln 


Vly = — & COS Py, ey = ep, | sin Bs | (3.9) 


i 
t 
; 
' 


: 
2 
: 


| 


: 
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f und 
f Bo = 2% —%. (3.10) 
Man geht deshalb am einfachsten so vor: 
(%y,. Vx.) > (€p,. %) mit UP 5, 
Bo nach (3.10) 
(,,> Pol > (% ty» ey) mit UP-4. 


Visierlinie Normflug 


Figur 1 


Zur Berechnung von ty, und e,,,, (Grundriss). 


Durch Division durch v, bzw. durch Absolutwertbildung erhalt man ty, bzw. 
€xy. Wegen (3.1) hat man dann: 


hy = hy (ty) = hy + 0, by - (3.11) 
-Damit sind alle Kenngréssen des Normfluges bestimmt. 


3.3 Berechnung der Eingangsfehler 


In kartesischen Koordinaten hat man als Eingangsfehler: 


Ax; = %y,— *y,» ANS Vin, In, 1 Aly hy Ty, » (3.12) 


wobei sich xy, , Yy, und hy, aus (3.1) fir t=] At berechnen lassen. Der Abstand 
der beiden Punkte M, und JN, ergibt sich dann zu 


Asq; = \ Ax? + Ay? + Ah? . (3.13) 


Eine dem vorliegenden Kgt. besser angepasste Darstellung der Horizontal- 
komponenten des Eingangsfehlers ist die Zerlegung in die Kursrichtung und 
die senkrecht dazu liegende Komponente. Die in Kursrichtung liegende Kom- 
ponente Aé ist dann in erster Linie fiir die Schwankungen der Horizontal- 
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kann man die Tatsache beniitzen, dass man sie durch Drehung des Koordina- — 
tensystems um den Winkel x aus Ax und Ay erhalt (Figur 2). Anstatt nun — 


Figur 2 


A& Eingangsfehlerkomponente in Kursrichtung (vorwarts positiv) ; 477 Eingangsfehlerkomponente 
senkrecht zur Kursrichtung (links positiv); 44 Héhenkomponente des Eingangsfehlers (nach oben 
positiv). 


aber die fiir diesen Fall iiblichen Transformationsformeln anzuwenden, ist es 
fiir die Programmierung einfacher, die Unterprogramme 5 und 4 anzuwenden: 


(Ax, Ay) > (4s,,~) mitUP5, 
p= px: 
(As,, p) > (4é, My) mit UP 4. 


Manchmal interessieren auch die Fehler der Eingangswerte «, A, e selbst: 
Ag = ty —aty, AA=Ay—dAy, Ae= ey — ey - (3.14) 


Dazu berechnet man die Kugelkoordinaten a, Ay, ey von N am besten mit 
Hilfe von UP 7 aus den kartesischen Koordinaten xy, yy, hy. 


4. Schussfehlerberechnung 
4.1 Approximation der Ballistik 


Die Ballistik wird beschrieben durch zwei Funktionen je zweier Variablen, 
zum Beispiel 


t(e, 4) = Geschossflugzeit,  e(e, 2) = Elevation. 


Diese beiden Funktionen sind entweder graphisch (Flugbahnkarte) oder tabel- 
lenmassig gegeben. Fiir digitale Rechenmaschinen kommen aber nur formel- 
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_ massige Ausdriicke fiir die beiden Funktionen t und ¢ in Frage, die nur elemen- 
tare Operationen enthalten. Sie wurden zu diesem Zweck durch Polynome 
approximiert. 

Als Approximationsbereich dient ein Rechteck in der (e, 4)-Ebene: 


(0 = ) emin = e = Cmax » Amin = A = higae (Ss 1600 A Io) i (4.1) 


Diesem Rechteck entspricht in der (e,, h)-Ebene ein Kreisringsektor, diirfte 
also dem Wirkungsbereich eines Flabgeschiitzes gut angepasst sein. Die lineare 
Transformation 

F=de—b, Ww=crA—da, (4. 2) 


in welcher a, b,c,d geeignete Konstanten sind, fiihrt das Approximations- 
rechteck (4.1) in das Einheitsquadrat itiber: 


ie aia Fe 1. (4.3) 


In diesen Einheitsvariablen € und 7 lauten die Polynome: 


n 


0 
= CC) Eel aa 
t= OO Eni, e= 


— nv 
1,4 


2, bY Ely). (4.4) 
15] 
Sie sind vom -ten Grade in jeder Variablen. 

Die Berechnung der Koeffizienten oy? und oi kann auf verschiedene Arten 
geschehen : 

Lassen sich t und ¢ an einigen Stiitzstellen geniigend genau bestimmen, so 
ist eine Tschebyscheff-Entwicklung vorzuziehen. Falls sie aber an vielen 
Stellen, aber weniger genau, bekannt sind, lasst sich eine Ausgleichsrechnung 
anwenden. Da t und « als Funktionen zweier Variablen stetig sind (man 
beschrankt sich ja bei der Flugabwehr auf einen eindeutigen Bereich innerhalb 
der Enveloppe der Geschossflugbahnen), lassen sie sich durch Ausdriicke der 
Form (4.4) bekanntlich beliebig genau approximieren. Die Erfahrung zeigte 
indessen, dass man mit 1 = 3 schon sehr brauchbare Approximationen finden 
kann. 

Eine den Formeln (4.4) entsprechende Approximation ist auch fiir die 
Derivation erforderlich. Wegen ihrer Kleinheit sind jedoch viel weniger Glieder 
nétig. Die Abhangigkeit vom Lagewinkel 4 kann gewohnlich vernachlassigt 
werden, so dass man nur noch eine Funktion von e allein hat, welche zum 
Beispiel folgende Form haben kann: 


b=e(d, ted), 


wobei d, und d, zu bestimmende Koeffizienten sind. 
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4.2 Berechnung des Treffpunktes 


Die graphischen Verfahren, derer man sich bei der Treffpunktberechnung ~ 
normalerweise bedient, sind einfach und iibersichtlich. Leider lassen sie sich 
aber bei digitalen Rechenmaschinen nicht anwenden. 

Man hat zunachst zwei Gleichungen. Die erste ist eine Vektorgleichung. 
Sie extrapoliert den Zielflug um die Zeit t + k in die Zukunft (k Ladeverzug): 


Yp =ty t+ (t+R) DV. (4.5) 


Diese Extrapolationsformel gilt natiirlich nur fiir gleichfoérmigen Flug, der uns 
hier allein interessiert. Bei ungleichf6rmigem Flug kann man sie entsprechend 
modifizieren, ohne dass die folgenden Uberlegungen ihre Giiltigkeit verlieren. 
Die zweite Gleichung liefert die Geschossflugzeit t als Funktion des Ortsvek- 
tors 77 des Treffpunktes T [siehe (4.4)]: 


hes Galle) ort (4.6) 


Durch Einsetzen der Zielfluggleichung (4.5) in die ballistische Gleichung (4.6) 
erhalt man eine einzige (nichtlineare) Gleichung fiir die Unbekannte T: 


tT=T(ry + (T+) 0). Cie 


Die Lésung kann zum Beispiel nach dem Newtonschen Iterationsverfahren 
erfolgen ([{3], § 27 C). Dieses besteht bekanntlich darin, dass man aus einer 
Naherung 1, fiir die Wurzel der Gleichung 


f(z) =0 (4.8) 


eine bessere Naherung 7;,, durch die Formel 


erhalt. In unserem Falle folgt durch Vergleich von (4.7) und (4.8): 
f(t) =t—t(ry + (t +h) 9d). (4.10) 
Mit Hilfe der Abkiirzung 
t(ry + (t; +k) ¥) = tF (4.11) 


hat man 
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_ Durch Differentiation von (4.10) folgt ferner (Anwendung der Kettenregel fiir 
_ eine Funktion mehrerer Variablen) : 


f(a) =1- (v- grad; t) , 


wobei der Ausdruck in runden Klammern das skalare Produkt der beiden 
_ Vektoren v und grad, t bedeutet. Das Einsetzen der beiden letzten Ausdriicke 
in die Formel (4.9) liefert als Iterationsformel : 


eae 4.12 
Tea = Ue = : ; 
"3 ie ies (v- grad, Tt) ( ) 


Als nullte Naherung t, wird man natiirlich die vom Kegt. gelieferte Tempierung 
' nehmen. Diese Naherung ist schon relativ gut, so dass die nachste Naherung tT, 
in fast allen Fallen schon geniigend genau sein wird (quadratische Konvergenz). 
Die Iteration wird abgebrochen, sobald das zweite Glied in (4.12) geniigend 
klein wird. Das skalare Produkt 


(v - grad t) (4.13) 


im Nenner liegt in folgenden Schranken, da der Vektor grad t in einer Vertikal- 
ebene durch das Kegt. liegt und ungefahr in Visierrichtung vom Kgt. wegweist : 


[ <0 im Anflugbereich, 


ja = (0 wenn 7 = W, 
(v - grad t) = 1 } ; (4.14) 
> 0 im Wegflugbereich, 


[ = 1 bei der Erreichbarkeitsgrenze des Wegfluges. | 


Die Erreichbarkeitsgrenze des Wegfluges liegt dort, wo bei der konventionellen 
Flugabwehr infolge Geschwindigkeitsverlust des Geschosses die letzte Treff- 
moglichkeit besteht. Aus (4.14) folgt, dass der Nenner von (4.12) praktisch 
immer positiv ist. 

Man kann das skalare Produkt (4.13), das ja durch die Ableitung von 
(4.7) nach t entstanden ist, durch einen Differenzenquotienten ersetzen: 


is Tt — 

(v- grad; t) & eee ; (4.15) 
wobei analog zur Abkiirzung (4.11) 

t (ty + (t; +8) 0) = 25" (4.16) 


gesetzt wurde. Man erhalt dann durch Einsetzen in (4.12) eine weitere Itera- 
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tionsformel, welche keine Ableitungen mehr enthalt: 


(«* a T;)? 


$2 GFT. 


Tin = Vi 


Diese Formel wird in der Mathematik nach STEFFENSEN benannt ((3], § 27,7) 
und dient ganz allgemein zur Konvergenzverbesserung von Iterationen [5]. 
Auf die aus (4.7) folgende «primitive» Iteration 


ee * 
Tin = T; 


sei hier nicht eingegangen, da sie in den meisten Fallen zu wenig rasch konver- 
giert [lineare Konvergenz mit dem Quotienten (4.13) ]. 

Ist t; durch Iteration bestimmt worden, so lasst sich der Treffpunkt T 
durch Einsetzen von tr in die Formel (4.5) berechnen. Damit lassen sich auch 
ér und a, bestimmen. Fiir die ballistische Koordinate «7 ist noch die Beriick- 
sichtigung der Derivation notig: 


ap =ap +o. (4.18) 


4.3 Berechnung des Sprengpunktes 


In mathematischer Formulierung ist die Sprengpunktberechnung nichts 


(4.17) 


anderes als die Transformation von ballistischen Koordinaten in kartesische — 


(Abschnitt 2.3). Es handelt sich also zunachst um die Umkehrung der Poly- 
nome (4.4): Gegeben sind t, und ¢,, gesucht sind é, und »,. Zur Berechnung 
dieser Umkehrung wurde wiederum das Verfahren von NEwTon herangezogen, 
diesmal jedoch formuliert fiir die Auflésung der folgenden zwei Gleichungen 


([3], § 27, (16)): 


t(é, 1) =) Uers e(é, 1) = €,- (4.19) 


Aus der 7-ten Naherung (&;, ;) ergibt sich die (¢ + 1)-te Naherung durch 


Gi = 6: +48, a =m+4n;, (4.20) 


wobei A; und An, Lésungen des folgenden linearen Gleichungssystems sind: 


OT 
ys 
0§ 


aa IS foams ee A= te Tes 


de | 
=| Abe + sal, -An; = & — &;. | 


(4.21) 
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- Dabei sind 
| OT 
os 


(Ot 
ine 07) 


Oe 
A MOS 


Og 


Ti» €;, ’ 
i On 


a 


_ die Werte der betreffenden Funktionen an der Stelle (&;, 7;). 
Man ist auch hier wiederum in der Lage, als nullte Naherung Werte &), No 
zu nehmen, die schon relativ nahe der Lésung sind, namlich das Resultat der 
vorhergehenden Treffpunktberechnung. Eine einmalige Iteration geniigt in 
den meisten Fallen. Die Iteration wird abgebrochen, sobald beide Fehler 
T, — T; und e, — e; dem Betrage nach unter einer gewissen Schranke liegen. 
Nachdem &, und , als Resultat der Iteration vorliegen, berechnet man 
daraus durch Umkehrung der Formeln (4.2) die Gréssen e, und /,. Diese ge- 
statten dann die Derivation zu bestimmen (UP 11), womit man dann die von 
der Derivation freie Seite berechnen kann: 


C=, 0- (4.22) 


Damit sind die Kugelkoordinaten «,, A, und e, des Sprengpunktes bekannt. Die 
_ kartesischen Koordinaten ergeben sich daraus nach (3.2) (UP 6). 


4.4 Berechnung der Schussfehler 


Die durch den Ladeverzug bedingte Komplikation bewirkt, dass das Ket. 
die zu einem Sprengpunkt gehérenden ballistischen Koordinaten nicht alle 
drei im gleichen Zeitmoment liefert. Vielmehr erscheinen « und ¢ um die Lade- 
verzugszeit k spater als das zugehdrige t. Deshalb muss #/A? eine ganze Zahl 

sein, namlich gleich der Anzahl der Bilder, um welche « und ¢ verschoben 
werden miissen, um mit t zusammen die ballistischen Koordinaten «,, ¢,, T, des 
Sprengpunktes zu geben. 
Die Differenzen 


Ac=%,—6p, Ae=&—ep, AT=%— ty (4.23) 


sind Kgt.-Fehler (Abschnitt 2.3). Die ballistischen (und auch die kartesischen) 
Koordinaten von T werden dabei von UP 14 geliefert. Zur Berechnung der 
Schussfehler braucht man ferner die kartesischen Koordinaten %,, y,, 2, des 
Sprengpunktes (UP 13). Fiir den theoretischen Schussfehler gilt dann (UP 15): 


TS = As, =) (x, — %7)* + (Vs — Yr)? + (hy — hy)? (4.24) 


Zur Bestimmung des wirklichen Schussfehlers braucht man noch den Gegner- 
punkt. Nach dem in Abschnitt 2.3 Gesagten gilt fiir ihn: 


¥o=ty + (tj +h) 0. (4.25) 


ZAMP IXa/24 
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Dann ergibt sich: 


GS = As, = (x, — %)? + (¥. — Ye)2 + (A, — Ng)? - (4.26) 


- 
| 


Bekanntlich ist im Anflugbereich der wirkliche Schussfehler grésser als der — 


theoretische, im Wegflugbereich kleiner; nur wenn der Treffpunkt in einer 
gewissen Umgebung des Wechselpunktes ist, kann beides vorkommen. 


5. Fehlerverteilung 


Die Berechnung der Fehlerverteilung ist relativ einfach. Sie wurde hier nur ~ 


durchgefiihrt fiir die Fehler As,, As, und As,,, welche ja nur positiv oder Null 
sein kénnen. Fiir jeden dieser drei Fehler werden 50 Zellen im Speicher der 
Rechenmaschine reserviert. Diese werden vor Beginn der Auswertung jedes 
Fluges mit Hilfe von UP 17 mit Nullen gefiillt. Ist dann im Laufe der Auswer- 
tung ein Fehler berechnet worden, so liefert seine Grésse die Adresse derjenigen 
Zelle, deren Inhalt um Eins erhéht werden soll. Am Schlusse hat man dann in 
jeder dieser Zellen die Anzahl derjenigen Fehler, deren Grésse der Adresse der 
betreffenden Zelle entspricht, das heisst, man hat die Dichteverteilung des 
Fehlers. 

Die Summenverteilung erhalt man daraus durch einfache Summenbildung. 
Diese sowie das Drucken der Summenverteilung wird von Hauptprogramm 3 
besorgt. Gedruckt wird folgende Funktion: 


A(As) ist die Anzahl derjenigen Fehler, welche kleiner sind als As +h 
(A Schrittweite), und zwar 


A,(As) fiir den Eingangsfehler As, , 
A,(As) fiir den theoretischen Schussfehler As, , 


A,,(AMs) fiir den wirklichen Schussfehler As,, . 


6. Zusammenstellung der verwendeten Rechenprogramme 
6.1 Unterprogramme 


Mit Ausnahme der Quadratwurzel, fiir welche das in der Rechenmaschine 
eingebaute UP beniitzt wurde, wurden alle UP selbst hergestellt. 

UP 1: Berechnung der Quadratwurzel. 

UP 2: Berechnung von sin x, auf 4 Stellen genau. x ist dabei in A °/5, gegeben 
und braucht nicht auf das Intervall 0 < x < 6400°/9 beschrankt zu 
sein, 

UP 3: Berechnung von arctg|x|, wobei x beliebig aber | «| +00. Das 
Resultat wird in A °/o9 geliefert mit 1/19 /o) Genauigkeit. 
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BUP «4; 


UIP Se 


UP 
SUP 


SIO 


BU P 9: 
ty PA0: 


PUP 11: 
UP 12: 


iP 13: 


UP 14: 


WPA: 


WP 16: 
WP 17: 
Wwe LS: 


Transformation von ebenen Polarkoordinaten in kartesische. cos x 
wird dabei als sin (« + 1600 °/o)) gerechnet. Beniitzt UP 2. 
Transformation von ebenen kartesischen in Polarkoordinaten, wobei 
der Winkel zwischen 0 und 6400 °/o, geliefert wird. Beniitzt UP 1 und 
Wheres 


: Transformation von Kugelkoordinaten in kartesische. Beniitzt UP 4. 
: Transformation von kartesischen in Kugelkoordinaten. Beniitzt UP 5. 
WP 28: 


Lineare Extrapolation in kartesischen Koordinaten um eine gegebene 
Zeit. 

Auflésung eines linearen Gleichungssystems zweiter Ordnung. 
Berechnung der ballistischen Polynome (4.4) mitsamt je den ersten 
partiellen Ableitungen (7 beliebig). 

Berechnung der Derivation. 

Transformation von kartesischen in ballistische Koordinaten. Beniitzt 
WP 710 und: 

Transformation von ballistischen in kartesische Koordinaten (Spreng- 
punktberechnung nach Abschnitt 4.3). Beniitzt UP 6, 9, 10 und 11. 
Treffpunktberechnung nach Abschnitt 4. 2. Liefert die kartesischen und 
die ballistischen Koordinaten des Treffpunktes. Beniitzt UP 4, 5, 8 
und 12. 

Berechnung des Abstandes zweier Punkte, deren kartesische Koordi- 
naten gegeben sind {zum Beispiel (4.26)]. Beniitzt UP 1. 

Verwertung eines Fehlers fiir die Dichteverteilung (Abschnitt 5). 
Nullsetzen der Dichteverteilung eines Fehlers. 

Trennung der drei in einem Wort zusammengefassten Koordinaten. 


Mit UP 18 hat es folgende Bewandtnis: Da ein Flug bis zu 200 Bilder haben 
kann, wurden zwecks Einsparung an Speicherzellen alle drei Koordinaten 
(ballistische oder Kugelkoordinaten) eines Punktes in einem 14stelligen Wert 
zusammengefasst und in einer einzigen Zelle gespeichert. Deshalb brauchten 
fiir die Eingangswerte und die Ausgangswerte nur je 200 Zellen reserviert zu 
werden. 

UP 2 und UP 3 (sin x und arctg x) verwenden Polynomapproximationen 
der betreffenden Funktionen in beschrankten Intervallen (siehe zum Bei- 


spiel [4]). 


EP ls 


6.2 Hauptprogramme 


Berechnung des Normfluges und Drucken von dessen Kenngréssen 
sowie Berechnung und Drucken der Eingangsfehler. Verlangt am An- 
fang die Eingabe (Tastatur) von Bildintervall At, Anfangsnummer 7, 
und Endnummer j,, der bei der Glattung zu beriicksichtigenden Bild- 


folge. 
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HP 2: Berechnung von Treff-, Spreng- und Gegnerpunkt, Drucken der Kgt.- ; 
und Schussfehler. Dabei wird der zuletzt mit HP 1 berechnete Norm- | 


flug verwendet. Verlangt die Eingabe des Ladeverzugs k. 


HP 3: Drucken der Summenverteilung von Eingangsfehler, theoretischem und ~ 


wirklichem Schussfehler. Die Dichteverteilung dieser Fehler wird 
bereits in HP 1 bzw. 2 berechnet. 

HP 4: Zur Kontrolle von UP 14. Verlangt Eingabe von «, A, e, %, v;, Vv, und k 
sowie eines Naherungswertes fiir 7. Liefert ¢,7, hy, tr, 6, x7, €r fiir den 
Fall k = 0 sowie tT, fiir das eingegebene k. 

HP 5: Zur Kontrolle von UP 13. Verlangt Eingabe von t; und és; sowie 
Naherungswerte fiir e,; und Ag. Liefert es, As, egs und hg. 

Zu HP 4 und HP 5 ist zu bemerken, dass vorausberechnete Ergebnisse 
(Einzelpunkte) schon langst in den Kontrolltabellen existieren, die zur statischen 
Priifung der Kgt. verwendet werden. Anhand dieser Tabellen kann deshalb die 
Treff- und Sprengpunktberechnung nachgepriift werden. 

Fiir das ganze Rechenprogramm sind total etwa 600 Befehlszellen erforder- 
lich, was iiber 1000 Befehlen entspricht. Insgesamt braucht es fiir Befehle, 
Konstante, Parameter und Variable im Maximalfall 1400 Zellen. Ferner wurden 
alle neun Indexregister verwendet. 


7. Schlussbemerkungen 


Die vorstehenden Darlegungen zeigen in groben Ziigen, wie die maschinelle 
Auswerterechnung bei der Flugabwehr durchgefiihrt werden kann. Leider ver- 
bietet es der vorhandene Platz, auf weitere Einzelheiten einzugehen. 

Man findet auch in diesem Falle wiederum die Tatsache bestatigt, dass 
programmegesteuerte elektronische Rechenmaschinen wesentlich andere Re- 
chenmethoden verlangen, als sie sonst iiblich sind. Hier betrifft es vor allem: 

Verwendung von Approximationen statt Tabellen; 

iterative statt graphische Auflésung von nichtlinearen Gleichungen. 

Ist das Rechenprogramm einmal gemacht (wie es hier zum Beispiel der 
Fall ist), so lassen sich damit beliebig viele Fliige auswerten. Die reine Rechen- 
zeit wird gegentiber den bisherigen Methoden schatzungsweise etwa um den 
Faktor 50 erniedrigt (Stunden statt Wochen). Dadurch wird es méglich, in 
kiirzerer Zeit wesentlich mehr statistisches Material zu gewinnen als bisher. 
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Résumé 


L’auteur montre dans les grandes lignes comment on procéde a l’analyse des 
épreuves optiques d’appareils de commande de tir (appareils directeurs) pour la 
_ défense terrestre contre avions sur une calculatrice digitale électronique dirigée 
par programme. I] donne entre autre des méthodes de calcul pour 
laplanissement des valeurs d’entrée, 
la représentation de la ballistique, 
le calcul du point d’impact, etc. 
Le programme de calcul pour l’analyse complete comprend dans le cas présent 
5 programmes principaux et 18 sous-programmes. 


(Eingegangen: 19. November 1957.) 
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A Method for Studying the Diffusion of Silver Iodide Particles 
in the Atmosphere by Means of I'*"1) 


By Rarer CarRreRAS-Patxot and RAYMUND SANGER, Ziirich?) 


1. Introduction 


Although the use of silver iodide vapours to stimulate either rain or various 
other meteorological phenomena is nowadays common practice, it runs up against 
an unknow factor that is crucial: where does the vapour go after leaving the 
generator ? 

A variety of authors have studied the question by employing a freezing chamber, 
which allows the icing nuclei to be counted on the basis of the ice crystals they 
have formed about themselves. In this way it has been possible to make clear that 
the number of nuclei is considerably increased in the neighbourhood of the generator 
and even over large distances in the direction of the wind. ; 

In the work presented here we have attempted to check the presence in the air 
of silver iodide, without relying upon those of its properties which might involve 
ambiguity. To this end we have tried marking the silver iodide with radioactive 
iodine 131. The originality of the method lies, we believe, in the use of relatively 


1) Paper presented at the 3rd International Symposium on Atmospheric Condensation Nuclei, 


Cambridge, England, July 16-18, 1958. 
2) Swiss Federal Institute of Technology, Department of Physics. 
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minute levels of activity, so that no danger attaches to working with the product © 
prior to its dispersal in the atmosphere. 
On the other hand this also means that we have to perfect measurement and © 
detection of extremely slight activity, of the order of 10-1’ curies per litre of air, 
which corresponds to one disintegration per litre per month! 
Let us follow, then, the particles from their production to the measurement of 
their activity: 


2. Production of the Vapour 


The process is well-known: the silver iodide marked with iodine 131, which we 
yet from Amersham, with an addition of sodium iodide is dissolved in acetone and — 
the mixture is burnt in a propane flame. Very fine dispersion is achieved in this 
way, the particles having a diameter of several hundredths of a micron (Figure 1). 


Figure 1 


Particles of silver iodide seen in the electronic microscope. 


3. Taking the Sample 


The extreme fineness of the particles makes very exacting demands on the 
filter which is to retain them. Above all it must offer sufficient flow for the slight 
differences of pressure that are available: the filter will, in fact, be mounted oan 
aircraft and it is only the dynamic pressure of the air which comes into play in 
bringing the air into the filter. The choice was decided in favour of a paper similar 
to that used in Swiss Army gas masks. 

This filter is placed in two airtight collecting tubes designed to be fitted under 
each wing of the aircraft (Figures 2 and 3), The volume of air filtered is measured 
by integrating the flow, which is recorded from the resistance in a platinum wire 


raised to red heat and cooled by the current of filtered air. Calibration is carried 
out in a wind tunnel, 


a. 
a 
—: 
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es From the way in which the filter is folded and then positioned in the collecting 
__tubes, it is possible to achieve a filter surface of nearly 2 m? in each collector with 
an aerodynamic section of only 2-5 cm?. 


ah ROM CEN Ba 


Figure 2 
Collecting tube. 


Figure 3 


The aircraft equiped with a filter. 


4. Concentration of the Sample 


After about 10 m® of air have been filtered, the traces of silver iodide must be 
extracted from the paper, and for this we first tried a biological method. 

It is an established fact that if iodine is given to mammals, a considerable 
proportion will become concentrated in the thyroid gland, which is generally tiny. 
We chose to use rabbits, even though these animals are known to be fussy about 
their food, Despite this psychological difficulty we managed to get them used to a 
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diet consisting of 50% filter paper and 50% cabbage leaves. Unfortunately, while - | 
the concentration in the thyroid was most satisfactory with regard to the sodium | 
iodide, the yield of silver iodide was extremely small, as a result probably of its — 
insolubility and its annoying tendency to form complexes of indeterminate structure. | 

It was necessary, therefore, to revert to a very much longer chemical process. 
In outline it is as follows: the filter paper, together with 0-2 mg of nonradioactive 
silver iodide, is calcinated slowly in the presence of sodium carbonate. The iodine 
is converted into sodium iodide which is then dissolved in water. After adjustment _ 
of the pH value around 7, silver nitrate is added to the solution, which gives a 
precipitate of silver iodide containing the activity of the iodine 131 from the filter. 

This precipitate is filtered, and the dry filter, with a surface area of about 10 cm?, 
is placed in a watch glass, into which is poured 2 cm? of acetone and 0-15 mg of 
sodium iodide. Thus we have the activity of the iodine from the filtered air concen- 


trated in these 2 cm? of solution. 


—_ 


5. Measurement of the Activity 


The radioactivity of this solution is so slight, however, that even with the best 
Geiger-Miller counters or scintillation counters one would not obtain more than 
several dozen impulses a day, whereas the background is around a dozen impulses a © 
minute. The method which has been used is based on the capacity of certain 
nuclear plates, such as Ilford G5, to register electron traces of some hundreds keV ; 
the beta electrons of iodine 131 are of this order of energy. 

The activity of the residue must, then, be concentrated into the smallest — 
possible surface, in order to achieve a maximum density of traces, which will | 
contrast most clearly with the inevitable background of the plate. 


Figure 4 


The measuring block with the active source in the middle. 


For this purpose the active solution is evaporated in a kind of crater hollowed 
out of an anticorodal block. The deposit is carefully detached from the walls and 
compressed at the bottom to form a small disc 1mm in diameter and 1/2 mm in 
height (Figure 4). This disc is then placed on a nuclear plate during two radioactive 
periods; for iodine 131 this is fifteen days. The plate is then developed. 
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| If we reckon that five electrons out of six will not reach the emulsion, and if we 

a have filtered a volume of 10 m# of air, then an activity of 10-17 c per litre of air will 
give about a thousand traces per square millimetre (Figure 5). 

However, if all possible precautions are taken (the newest plates are used, 

_ contamination is suppressed, etc.) the number of traces from sources other than the 


Figure 5 


Traces corresponding to total activity of 20, 1, 1/50 pe. 


iodine 131, and which might be confused with those we are looking for, amounts to 
several dozen per square millimetre. This number is relatively small because, if we 
make a visual count (and this has proved the only easy method), then parasite 
_ traces due to alpha rays, such as those caused by high energy particles coming 
from cosmic rays, can be clearly recognized and disregarded. 


Conclusions 


This method has only recently been perfected and we have not yet been able 
to conduct many experiments. The first results obtained with a very weak emission 
of 0-1 uc/s have enabled us to detect silver iodide at a distance of several hundred 
metres from the generator. With an emission one hundred times stronger, of the 
type that is normally planned, it should be possible to detect silver iodide up to a 
dozen kilometres away, by exploiting the maximum sensitivity of our method; 
this would be its furthest range. 

By reducing the dimensions of the active surface, sensitivity could be increased 
by a factor of 100 without great difficulty ; but then the iodine which arises from the 
nuclear explosions in the world, and which is present most irregularly, will distort 
the results. 

Although our point of departure is quite different from that of authors who use 
a freezing chamber for the same research, we finally encounter the same difficulty : 
is what we are detecting the silver iodide which we produced on the ground, or is 


it something else ? 
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| 
If one is content to utilize such small amounts of activity that there is no © 
practical problem of protection from radiation, then radioactive marking of silver 
iodide to study its diffusion is a process which, with all its advantages, snags and 
natural limitations, can be added to existing methods, even though it has not 
produced the revolutionary results that we are perhaps overinclined to expect in | 
our present fascination with the wonders achieved by radioactive isotopes in other 
fields. : 
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Zusammenfassung 


Es wird eine neue Methode entwickelt zur Erfassung der Diffusion des zur 
kinstlichen Wetterbeeinflussung verwendeten AglI-Rauches in der freien Atmo- 
sphare. Der durch Verwendung von I! radioaktiv gemachte AgI-Rauch wird in 
speziellen Papierfiltern, die am Messflugzeug befestigt sind, aufgefangen; mittels 
eines chemischen Prozesses wird das AgI aus den Filtern gelést und sodann auf 
eine Kernphotoplatte gebracht, wo die Zahl der B-Spuren ein Mass fiir die Anzahl 
der in der Luft gefundenen AglI-Partikel darstellt. 


(Received: May 12, 1958.) 


The Solution of Cubic Equations by Iteration 


By WILLIAM ANDREW CopPpEL, London, England?) 


1. A very simple iterative method of factorizing polynomials has been suggested 
by Lin [1]?). Applied to the determination of a linear factor, that is, to the determi- 
nation of a root of the polynomial equation p(*) = 0 it takes the form of calcu- 
lating a sequence of successive approximations (x, according to the rule 


poems 
ea = TBO) = Pl) 


1) Department of Mathematics, Birkbeck College. 
*) Numbers in brackets refer to References, page 383. 
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_ Since {p(0) — p(0)}/% is a polynomial of degree one less than p(#) itself the number 


_ of operations needed to produce a new approximation is reduced to a minimum. 


On account of its simplicity the method can be less time-consuming than other 
methods, such as NEwron’s, with a nominally higher order of convergence. The 
disadvantage of the method lies in its uncertain convergence. The only criteria 
which have been proposed (see, e.g. [2]) establish convergence under certain 
conditions provided the starting-point lies in a sufficiently small neighbourhood of 


_ the root. But as the root has yet to be found the interval in which the starting-point 


may be taken is initially unknown. A criterion which does not suffer from this 
defect is obtained in the following for the simplest case of a real cubic equation. 
It turns out that convergence is closely connected with the position of the real 
parts of the roots, and in the type of problem for which the method was originally 
proposed this is known in advance. 


2. Let 
p(*) = + a27°+be%+c=0 
be the given cubic equation. If we put 


c 


NB core: emer ean 


_ then the successive approximations to a root of p(x) are determined by the recur- 


rence relation %,,, = [(%,). 

The function /(v) is defined for all real values of w if and only if a? < 4b. We 
will suppose this condition satisfied. If c = 0 then f/(7) = 0 and the iteration process 
converges immediately to the root ¥ = 0 of p(w). We may therefore suppose c + 0. 
We can actually assume c < 0. Otherwise, if c > 0, we have only to take — ¥ as 
the unknown in order to obtain an equation of the same form with negative 
constant term. 

Under these assumptions /(#) is a positive, continuous function which increases 
from zero at x = —co to a maximum value at x = —a/2, and then decreases to 
zero again at = +00. We assume finally that f(x) < 7 at x = —a/2. That is, we 
assume $(— a/2) > 0 or, in terms of the coefficients, 


ti 1 
3 
a peal alicia 


This last condition ensures that all real roots of the equation /(¥) = x, which is 
equivalent to the given cubic equation, fall inside the interval (—oo, —a/2) in 
which f(x) is an increasing function. We can therefore apply to this interval the 
following elementary result (cf. [3]): 

Let f(x) be a continuous, increasing function in an interval « Sx SB in which 
a <f(x) <f. Starting from any point x» im (a, B) let (x,) be the sequence defined 
vecurvently by setting %»,, = f(¥,). Then the equation f(x) = x has a root situated on 
the same side of x) as x, and the sequence (%,) converges monotonically to the nearest 
such root. Y 

We deduce that every iteration sequence whose starting-point lies in the interval 
(—oo, —a/2) converges monotonically to a root of p(x). If p(x) has three real roots 
the sequence converges to the greatest or least root according as the starting value 
wv, is greater or less than the middle root. The restriction on the starting point 1s 
of no consequence, since +, = /(¥9) lies in the interval (0, —a/2) for arbitrary %». 
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It is remarkable that the three conditions which we have had to impose on the 
coefficients of p(w), namely 


} 
b> La, ol ab feeeah ands to0e (1) 


are equivalent to the single condition that all roots of p(x) have real part less than — 


—a/2. In fact all roots of p(v) have negative real part if and only if b > 0, ¢ > 0, 
and ab > [4]. For the cubic 


il 1 ; i. se) ol. 
p (x 5 a) =e ae + (8 pvlrte—zabt =a 


these conditions are exactly the inequalities (1). Since —a is the sum of the real 
parts of all three roots an equivalent formulation is that the real part of any root 
is less than the sum of the veal parts of the remaining roots. If we denote the three 
real parts by 71, %2, 73 (¥%1 S % S73) it follows that 


1 
ee 8 ip Sa Se Oe areal: 


i 


To make our iteration process completely automatic we can prescribe the 


initial point v7). For *, = —a/2 we will always arrive at the greatest real root of 
p(x), and for 0 < %,< —a/4 at the least real root. Since the roots may le any- 
where between 0 and —a/2 the best universal starting-point is 7) = —a/4. 


The inequalities (1) can be somewhat simplified. In fact since the last two 
inequalities imply a (b — a?/4) < 0 the first inequality can be replaced by a < 0. 
Thus (1) is equivalent to 


a<0, p(—j a) >0, and p(0) <0. (2) 


If p(x) has only one real root the first inequality actually follows from the other two. 
3. Consider now an arbitrary cubic p(x) and let 
Pls) =P (# — @) = #8 4+ (a— 3.d) x? 4+.---, 


The roots of p() are obtained by adding d to the roots of p(v). Therefore p(x) will 
be a cubic of the type just considered if all roots of p(*) have real part less than 
(d — a)/2. Thus the iteration process can always be applied with certainty when 
an upper bound for the real parts of the roots is known. 

If nothing is known about the position of the roots we can try to choose d so 
that p(x) satisfies the conditions (2). That is, we require 


1 (d—a 
toh ss mae (=>) >0, and p(—a) <0. 
The first condition can be omitted if it is known that not all roots of p(«) are real. 
Evidently there exists a value d, such that these requirements are met by alld = dp. 
Hence the iteration process can be applied with certain success to an arbitrary cubic 
equation. It is not difficult to show that we can always take 


dy =a+2+2max(la|, |d], cl), 


but when a, b, c have numerical values it is usually possible, and preferable, to find 
much smaller values for d. 
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It is perhaps worth remarking that with given numerical values for the coeffi- 
cients of p(x) the coefficients of (7) are in general most simply obtained by using 
_ Horner’s well-known scheme for ‘synthetic division’. 

There is one frequently occurring type of problem for which the present method 

seems especially appropriate. In the solution of linear differential equations with 
constant coefficients it is often known a priori, or is determined by a preliminary 
investigation, that all roots of the characteristic equation have negative real parts. 

_ Thus, in the case of a third order system, we can take d = a in the analysis above. 
With this value of d we have 


Ss P(x) = p (*¥ — a) = #8 — 2ax? + (224+ bd) 4+ c-—ab. 


The coefficients here are particularly simple and the formation of f(x) determines 
of itself if the stability conditions are satisfied: (7) must have positive coefficients 
and (7) must have negative constant coefficient. 
The equations originally encountered by Lin were of just this type and arose 
in a study of the effects of gusts on aeroplanes. Lin observed that in some cases 
where his method failed convergence could be brought about by a suitable change 
_ of variable v = x’ — d. The results which have been established here explain why 
_ this is so and give a practical prescription for the choice of d: take d = a. 
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Riassunto 


L’autore tratta la convergenza del metodo iterativo di Lrn nel caso di un’equa- 
zione cubica. Stabilisce una relazione fra la convergenza e la grandezza delle parti 
_ reali delle radici, e dimostra che una radice puo essere sempre trovata in questo modo, 
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Low Prandtl-Number Free Convection 


By Joun L. Grece and Epurar M. Sparrow, Cleveland, Ohio, USA.*) 


Introduction 


PRANDTL’s boundary layer concept was first applied to the analysis of free 
convection by Scumrpt and BECKMANN [1]?) in 1930. They considered the problem 
of an isothermal vertical plate situated in a large body of quiescent fluid whose 


1) NACA, Lewis Flight Propulsion Laboratory. . 
2) Numbers in brackets refer to Reference, pages 387. 
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temperature is different from that of the plate. The mathematical model and — 


coordinate system are shown in the sketch below. 

These first investigators focused their attention on air and obtained a series 
solution of the laminar boundary layer (ordinary) differential equations*) for a 
Prandtl number of 0-733. Later, numerical solutions of these same boundary layer 
equations were carried out for Prandtl numbers between 0-72 and 1000 by ScuuH [2] 
and OsTRACcH [3]. 

The increasing interest in liquid metals as heat transfer media has directed 
attention toward an altogether different range of Prandtl numbers. Technically 
important liquid metals have Prandtl numbers lying between 0-003 and 0-03. For 
this range, isolated numerical solutions of the laminar boundary layer equations for 


Gravity 
field 


Figure 1 
Coordinate system. 


free convection on an isothermal vertical plate have been contributed by OstTRacH[3] 
for Py = 0-01 and SuGawara and Micutyosui [4] for Py = 0-03. The aim of the 
present note is to summarize new numerical solutions which extend over the entire 
liquid metal range. Comparisons will be made with the previous work noted above 
as well as with an approximate Karman-Pohlhausen formulation due to EcKERT [5] 
and SgurtrRE [6]. We will deal here with the most important quantity of technical 
interest, namely, the heat transfer. Details of the temperature distributions must be 
omitted because of space limitations. 


Brief Review of Theory 


Using the mathematical ideas of boundary layer theory, Scumipt and BEcK- 
MANN introduced a new independent variable 7 by the transformation 


49? (1) 


= 1/4 
nC Lee Cm [fF To Ta) 
x 
where g is the acceleration of gravity, 6 the thermal expansion coefficient, and 
v the kinematic viscosity. Then, the stream function y and the temperature T were 
related in the following way to new dependent variables F and @ which are functions 
only of ». 


year Cae Fin) pT = OG T Tal Seles. (2) 


In terms of the new variables, the equations expressing conservation of mass, momen- 


3) See equations (3) and (4). 


H 
4 


teal 


dq 
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tum, and energy for laminar free convection on an isothermal vertical plate become 
PEP SEE — 2(F Ye + 9 = 0 (3) 


and 
/ 6” + 3(Pr) F 0’=0, (4) 


where the primes represent differentiation with respect to , and Py denotes the 
- Prandtl number. 
The physical boundary conditions require that both velocity components vanish 
_at the plate surface (vy = 0), and that the vertical velocity component vanish far 
from the plate as y > oo. It is also required that the temperature T take on the 
value T,, at the wall and T,, in the ambient fluid. The expression of these conditions 
in terms of the transformed variables is 


F=0 ; 

aed = 0 S > oo 

“ a 5 r Fe 5 
fe AESE 7 (5) 

6= 1 


The conditions for 1 + oo are always troublesome in any numerical treatment and 
can only be satisfied approximately. Inasmuch as the new solutions were obtained 
with a relatively fast electronic computer, it was possible to satisfy the boundary 
conditions for large 7 more closely than in previous work which was done on desk 
machines or on slower computers. This circumstance, coupled with the fact that 
smaller steps were taken in the numerical process, suggests that the present results 
may be more accurate than those obtained in the past. 

Numerical solutions of equations (3) and (4) subject to conditions (5) have been 
carried out for Prandtl numbers of 0:03, 0-02, 0-008, and 0-003, and the results will 
be used in the heat transfer computation below. 


Heat Transfer Results 


The heat transfer results are reported in terms of the local heat-transfer coef- 
ficient, h, and local Nusselt number, Nu,. Also involved are the Grashof number, 
Gr,, and the Prandtl number. These several quantities are defined as follows 


oe: q = h% — & B (Zw — Tyo) ae y= Cp M 
h — Pee Nu, = 3 GH = y2 , Pr= k > 


where q is the local heat flux. In terms of the variables of the analysis, a dimension- 
less representation of the local heat-transfer results may be written as 
Nu, (—0'(0)] = 


(Gr, Pr?)t/4 = (2 Py)il2 * 


The quantity 6’(0) is an abbreviation for [d0/dy],,9 and is found from the solutions 
of equations (3) and (4). 

A tabulation of heat transfer results based on existing numerical solutions of 
equations (3) and (4) for the low Prandtl number range is given in the table. It is 
immediately seen that the quantity Nw,/(Gr, Pr*)'/* is relatively independent of 
the Prandtl number, and this is the reason for its selection as a correlating group. 


fi 
5 
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i 
It is noteworthy that dimensional analysis predicts such a finding provided that | 
viscous forces are omitted from the momentum equation. 


Pr | Nup|(Gry Pr®)H4 Investigation 
0-03 0:5497 Present analysis 
0-03 0-555 See [4] 
0-02 0-5582 Present analysis. 
0-01 0-574 See [3] 
0-008 0-5729 Present analysis 
| 0-003 0-5827 * | | Present analysis 


To facilitate comparison with other solutions, the results of the present study 
have been plotted on Figure 1. Also shown are OsTRACcH’s. point at Py = 0-01 and 
SUGAWARA-MICHTYOSHI’s point at Py = 0-03. These latter results fall sightly above 
our results (within 1%); such deviations are to be expected from their earlier, less 
accurate solutions. The heat transfer result from the Karman-Pohlhausen formu- 
lation as carried out by ECKERT and SQUIRE is 


Niepineret 0-508 (7) 
(Gr, Pr?)t4 ~ (Pr + 0-952)14 * 


Equation (7) is also plotted on Figure 1 and lies from 7% to 12% below the present. 
numerical solution, the greater deviations occuring at the lower Prandtl numbers. 


ke — 
1 Ei ie 
0-6 Present analysis,aumerical solution | | 
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kéiméin-Pohthausen solution|58} 
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t+ — rf 


0-4 ai =I — 
0:003 0-004 0-006 0:008 0-01 0-02 0:03 0-04 
Prandtl number, Cp /t/k —— 
Figure 2 
Local heat transfer results for low Prandtl number fluids (hay = 4 hyp/8). 


\J [8], numerical solution; O [4], numerical solution. 


This agreement is remarkably good when one considers the relatively broad as- 


ace employed in carrying out the Karmén-Pohlhausen procedure for this 
problem, 


A 
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For engineering purposes, a simple representation of our results for the liquid 
- metal range is 
Nt 402 505 (Gye Py2\ts (8) 


This formula represents the numerical solutions to an accuracy of better than 3%. 
' For an isothermal plate of length L, it is easy to show that the average heat- 
transfer coefficient h,,, is equal to 4/3 of the local coefficient at x = L. 
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Zusammenfassung 


Eine Untersuchung der freien Konvektion wurde an einer senkrechten isother- 
men Platte bei Prandtlscher Zahl im Bereiche von 0,003-—0,03 (das heisst fliissige 
Metalle) durchgefiihrt. Dabei wurden zahlenmassige Lésungen der laminaren 

_Grenzschichtgleichungen erhalten. Fiir die 6rtlichen Warmeiibergange bei diesem 
Bereiche der Prandtlschen Zahl lassen sich die Resultate als Nu, = 0,565 (Gr, Pr?)+/4 
_ darstellen. 
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A Note on Magneto-Hydrodynamics of a Finite Rotating Disk 


By SHyAMAL Kumar Majumpar, Kharagpur, India’) 


In a recent paper STEWARTSON?) has considered the problem of steady rotation 
of a finite disk in an electrically conducting liquid in the presence of a magnetic 
field. In the solution of the problem, the radial component H,, of the magnetic field 
has been neglected in comparison with the ‘strong’ external field H imposed along 
the axis of rotation. In this note it is shown from considerations of geometric 
dimensions that for a rotationally symmetric steady hydromagnetic field the 
radial component of the magnetic field can always be neglected in comparison 
with the field H,, along the axis of rotation if the depth of liquid column d is small 
compared with the radius a of the disk. It is also shown that in LEHNERT’s?) 


experiments this condition is satisfied. 
1) Indian Institute of Technology. 


2) K. SrewarTson, Quart. Mech. appl. Math. 10 (1957). 
3) B. Leunert, Proc. Roy. Soc. Lond. [A] 233 (1955). 
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Proof. Since the disk is in steady rotation the relevant equations are 


curlE=0, (1) § 
cunlH=4a7, ; (2) 
J=ClE+ pox A)], (3}4 
Za 


uy 


seen 


Trough containing liquid the base of which is rotating about the vertical axis (OZ). a radius of the 
rotating disk; d depth of the liquid column. 


where the vectors and the constants (uw, o) have their usual meaning. From (2) 
and (3) 
culH=4a2o0[E+u(v x A)). (4) 


In a system of cylindrical polar co-ordinates (7, #, z) on account of rotational 
symmetry all dependent variables are independent of 6. Hence the radial component 
of (1) and the azimuthal component of (4) reduce to 


OF, _ 
Soe= a ? (5) 
0H, Onl yee 
From (5) and (6) one obtains 
OM 2 08H. 
“et Or Oz ”) 


Since a and d are representative lengths in the radial and axial directions respec- 
tively, equation (7) suggests the following relation between the orders of magnitude 


H, 4H, 
@ ad (8) 
or 
d /) 
H, = — H,. (8”) 


Thus if dja <1, H, is small compared to H,. 
In the experiment described by LEHNERT the mean diameter of the rotating 


copper rings is 7 cm and the depth of mercury is 6 mm; consequently the neglect 
of radial magnetic field is justified. 


I thank Prof. B. R. Seru for his kind interest in the work. 
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Zusammenfassung 


Es wird gezeigt, dass in elektrisch leitender Fliissigkeit bei zylindrischer Sym- 
_metrie die Radialkomponente des Magnetfeldes klein ist im Vergleich zur Axial- 

komponente, sobald die axiale Ausdehnung des Gefasses klein ist im Vergleich zu 
seinem Durchmesser. 


(Received: May 16, 1958.) 


Zur Frage der Charakterisierung stationarer Bewegungen 
in der Hydrodynamik 


Von HENNING MULLER, Mainz, Deutschland!) 


Nach einem Satz von HELMHOLTz und KortEwEc [1]?) ist die unter Einwirkung 
von konstanten Kraften eindeutigen Potentials zustande kommende Bewegung 
einer reibenden Fliissigkeit durch die Eigenschaft ausgezeichnet, dass die Energie- 
dissipation unter der Voraussetzung vernachlassigbarer Tragheitsglieder fiir jedes 
feste Gebiet im stationdren Zustand kleiner ist als im Fall irgendeiner anderen 
Bewegung mit denselben Randwerten der Geschwindigkeit. Fiir diesen Satz wird 
eine Verallgemeinerung gegeben und eine Anwendung diskutiert. 

Betrachtet werde die Str6mung einer inkompressiblen, reibenden Flissigkeit. 
Der Ausdruck fiir die Energiedissipation in der Zeiteinheit hat die Form 


OU; 
fo dt = [ris vg dt 


(wobei— wie auch im folgenden — in tiblicher Weise tiber die zweimal vorkommenden 

Indizes — hier 7 und k — zu summieren ist). 7;;, bedeutet den Reibungstensor, v, den 

Geschwindigkeitsvektor und wv, den Ortsvektor. Vergleicht man zwei Bewegungs- 

zustande v, und v,;— wobei v;, — vj; am Rande des Integrationsbereiches verschwin- 
den soll —, so kann man schreiben 


[ews dx — [ 2) dr -|® (vi, — v,) de + [2 dx (1) 


- oder, da ® unabhangig von seinem Argument, wegen 
Yin = 2 deg 


(b;, bedeutet den symmetrischen Teil des Tensors der Verschiebungsgeschwindig- 
keit grad; y und 7 den Koeffizientén der inneren Reibung) stets positiv ist, 


[een dt — [ %@) ie > [2 dx. 2) 


Die Funktion Q ist jetzt folgendermassen zu ermitteln [2] (der Tensor 7;;, gehdrt 
zum Stromungsfeld v;,) : 


(Tiaey ares OU; ; O(0;, — Ux) P 
ty Qdt= / (“ix ae: Vike ov, ("1% — Yer) — 5 wae | dt ; 
a t 


1) Institut fiir theoretische Physik der Universitat Mainz. ' 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 392. 
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wegen 
dv; ,. Uy 
"ik fe : a 
ist 
ov; OUR 
Vik My EE ee 
t 
und damit 


OU, ra) 
[oae=2frn (Qe — ae) a 
a i 


Schliesslich folgt wegen des Verschwindens von v; — v;,am Rande des Integrations- 
bereiches durch partielle Integration 


[oa=2fom—up = dt 
v 


Ws 3 OU; 07v;, 


v 
Ox; 04, O%; oh ov? 


Wegen 


und wegen der Voraussetzung der Inkompressibilitat ist damit 
[aar=2nf mv) Av,, dt . (3) 


Der Differentialausdruck Av, kann nun aus der Bewegungsgleichung entnom- 
men werden; wird einschrankend vorausgesetzt, dass v, einen stationaren Stré- 
mungszustand kennzeichnet, so lautet sie 


OU; 
Co ee Ad eee (4) 
at} ‘k 


wobei ¢ natiirlich das Potential der ausseren Krafte, g die Dichte und # den Druck 
bedeuten. (3) in Verbindung mit (4) ergibt fiir 2 jetzt folgendes: 


Wegen der Ubereinstimmung von v, und v; am Rande o des Integrationsbereiches 
ist 


und natiirlich 


also 


oder 
: ra) 
[een den? e | (oy, — vf) 0, 2% de > [ dz. (5) 
tl 


Diese Extremalaussage — im Fall 9 = 0 besagt sie, dass in irgendeinem festen 
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Gebiet unter den gegebenen Voraussetzungen der stationdren Bewegung einer 
_ Fliissigkeit ein Kleinstwert der Energiedissipation entspricht — soll jetzt auf eine 
spezielle Bewegung angewandt werden, und zwar auf das an anderer Stelle unter- 
suchte Verhalten eines in einfacher Weise schematisierten Makromolekiilmodells in 
einem inhomogenen Strémungsfeld [3, 4]. Dass man mit Hilfe von (5) hierbei iiber- 
‘haupt zu einer einerseits richtigen und andererseits nichttrivialen Aussage kommen 
kann, ist keineswegs selbstverstandlich. Wahrend (5) sich auf eine reine Fliissig- 

_ keitsstr6mung bezieht, wobei von irgendwelchen, in gewissem Sinne mitschwim- 
menden Partikeln gar nicht die Rede ist, basieren die Uberlegungen, die sich an 
(6) anschliessen werden, auf der unter dem Einfluss Stokesscher Reibungskrafte 
zustande kommenden Bewegung eines in eine inhomogen, aber ungestért stromende, 
reibende Fliissigkeit suspendierten, in geeigneter Weise schematisierten Ké6rpers. 

Ausgangspunkt ist die Differentialgleichung 


@8—12275s a? (qsin?d— 8) =0. (6) 


Ihr liegt die Vorstellung der zweidimensionalen Bewegung einer Fliissigkeit (Zahig- 
keit 7), charakterisiert durch einen Geschwindigkeitsvektor v,, = (0, ¥q) (0 und 
xq bedeuten die Komponenten des Geschwindigkeitsvektors in bezug auf ein 
kartesisches Koordinatensystem mit den Achsen ¥ und y, q ist eine Konstante) und 
eines Hantelmodells (Tragheitsmoment @) zugrunde, das aus zwei schweren 
-Kugeln (Radien s) besteht, deren Schwerpunkte den festen Abstand 2a haben. 
# kennzeichnet die Lage der Hantel in bezug auf die y-Achse. Krafte werden durch 
Stokessche Reibung tibertragen. Die Differentialgleichung (6) beschreibt mithin 
_ die Bewegung des Hantelmodells, die dieses unter dem Einfluss des durch v, , gege- 
benen inhomogenen Strémungsfeldes ausfiihrt. - Durch eine Transformation 


d= q- 2(8) 
und mit der Abkiirzung 
O” 
¥ 122 sa? 
geht (6) tiber in 
yee —sin?e+z2=0. (7) 


Unter der Annahme einer periodischen Lésung 
2(B) = z (8+ a) = Zpep(0) 


— die Diskussion Von (7) etwa mit Hilfe des Isoklinenverfahrens bestatigt die aus 
physikalischen Griinden naheliegende Vermutung, dass sich alle Bewegungen des 
Modells schliesslich einer periodischen Bewegung asymptotisch nahern — liefert die 
Integration von @ bis # + x (tritt # als Integrationsvariable auf, so wird an seiner 
Stelle der Buchstabe ¢ geschrieben) 


Diese Aussage bezieht sich also, um es zu wiederholen, auf eine periodische Bewe- 
gung des Modells im ungestorten Strémungsfeld, auf eine Bewegung, der «im 
Grossen» der Zustand der Stationaritat entspricht. - Nun der Ausdruck fiir die 


“tly 


a 
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Energiedissipation: Bedeutet u die vektorielle Stromungsgeschwindigkeit am Ort — 
des Schwerpunktes einer der beiden Kugeln, die — starr miteinander verbunden — — 
das Teilchenmodell darstellen, und bezogen auf diesen Schwerpunkt, so ist die 
in Reibungswarme umgesetzte Energie 


WB 12ams [ude = 121s | ue dt 
mit 

u2 = a? (qg? sin?d + # — 2q bsin?d) , 
und zwar natiirlich sowohl fiir den Fall stationarer als auch fiir den Fall nicht- 
stationdrer Bewegung. Es liegt jetzt nahe, die in (5) aufgetretene Energiedissipa- 


tion mit der durch die Bewegung des Hantelkérpers in der Zeiteinheit verzehrten ~ 
Energie zu identifizieren, das heisst zu schreiben 


© = 12275 a7 (gq? sin?6 + & — 2q Gsin?d) . 


Die Ungleichung (5) lautet damit, natiirlich unter Verzicht auf die raumliche 
Integration, folgendermassen: 


676.98: (2% 22 pg = 2 Si A ie) io eee) pee nor ek oe 


Division durch z — z,,, und Integration nach @ tiber einen Bereich a liefert 


O+7 


[e+ Aon) aS x, 
—d 


wenn z— 2,,,S$ 0, und diese Ungleichung wird erfiillt, das heisst, sie stimmt fiir 
alle méglichen durch die Differentialgleichung (7) erfassten Bewegungen des be- 
trachteten Modells im Rahmen der angegebenen Schematisierung. 
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Summary 


HELMHOLTZ and KorrEewec propose that the steady motion of a viscous fluid 
under constant extraneous forces having a single-valued potential dissipates—for 
any given region and assuming that inertia terms in the dynamic equations can be 
neglected—less energy than any other motion with the same values of velocity at 
the boundary.—A generalization of this proposition is here given, and an appli- 
cation discussed. The application deals with the motion of a simple macromolecule 


model in an inhomogeneous field of flow—a motion caused only by the influence of 
SToKEs’ friction. 


(Eingegangen: 14, Marz 1958.) 
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Handbuch der K4ltetechnik. Herausgegeben von R. Pranx. Achter Band: 
H. Hausen, Evzeugung sehr tiefer Temperaturen. Gasverfliissigung und Zerlegung von 
Gasgemischen (Springer-Verlag, Berlin 1957). 412 S., 307 Abb.; DM 72.-. 

Das zwolfbandige Werk R. PLanxs, dessen achter Band jetzt erschienen ist, 


_ bedarf wohl keiner nahern Beschreibung, da es dem Physiker und Ingénieur, der 
_ sich mit der Kaltetechnik befasst, wohlbekannt ist. 


Die Erzeugung sehr tiefer Temperaturen, womit sich H. HausrEN im vor- 


_liegenden Band eingehend auseinandersetzt, hat in der letzten Zeit auf allen 


Gebieten der Technik sehr an Bedeutung zugenommen. 
Der erste Teil befasst sich mit der Theorie der Gasverfliissigung und der Zer- 
legung der Gasgemische. Die dazu n6étigen thermodynamischen Beziehungen 


sind in knapper und iibersichtlicher Form behandelt, gefolgt von Zustandsdia- 


grammen der tiefsiedenden Gase. In angenehmer Weise wird der theoretische 
Teil durch kurze Beschreibungen typischer Anwendungen der Technik aufge- 
lockert. Ferner ist der Theorie der Warmetibertragung im Gegenstrom, die fiir 


die Ausfiihrung solcher Apparate und Anlagen von grundsatzlicher Bedeutung 


ist, ein separates Kapitel gewidmet. 

Der zweite Teil befasst sich mit der Ausfiihrung und dem Betrieb der Ver- 
fliissigungs- und Zerlegungsapparate sowie mit der Anwendung der gewonnenen 
Fliissigkeiten und Gase. Dabei streift der Verfasser die in der Technik bewahrten 
Ausfiihrungen von Kompressoren, Expansionsmaschinen, Warmeaustauschern 
und Rektifikationskolonnen. Uber den Betrieb solcher Anlagen liegen verschie- 


_ dene praktische Hinweise, Betriebsergebnisse und Kostenrechnungen vor. 


Abschliessend kann vorliegender Band, als ein wirklich brauchbares und um- 
fassendes Buch der Tieftemperaturtechnik, empfohlen werden. Nicht vergessen 
sei, auf die zahlreichen Literaturangaben hinzuweisen. H. GRUTER 


Handbuch der Laplace-Transformation. Von Gustav Doetscu. Bd. 2: 
Anwendungen der Laplace-Tvansformation, 1. Abteilung, 436 S., 48 Fig.; Ganz- 
leinen Fr./DM 56.15, broschiert 52.— (1955). Bd. 3: Anwendungen der Laplace- 
Transformation, 2. Abteilung, 300 S., 23 Fig.; Ganzleinen Fr./DM 40.-, broschiert 
36.— (Birkhauser Verlag, Basel und Stuttgart 1956). 

Unter dem Titel der beiden Bande kénnte man zunachst ein Werk vermuten, 
welches die Lésung von technischen Problemen mit Hilfe der Laplace-Transfor- 
mation behandelt, nachdem in Band 1 die theoretischen Grundlagen bearbeitet 
wurden. Indessen ist hier der Begriff der Anwendung in einem sehr viel weiteren 
Sinne zu verstehen; es werden namlich Anwendungen der Laplace-Transforma- 
tion auf Probleme aus den verschiedensten Gebieten der Mathematik behandelt. 
So gibt das ausserordentlich reichhaltige — stellenweise fast zu reichhaltige — 
Material nicht nur eine Ubersicht iiber die Anwendungsméglichkeiten der Laplace- 
Transformation, sondern ist auch eine Fundgrube fur funktionentheoretische 


Beziehungen. ; ‘ 
Der Inhalt ist wie folgt in 7 Teile gegliedert (Teil 1-3 in Band 2, Teil 4-7 in 


Band 3): 

Der erste Teil behandelt die Vielfalt der Anwendungsméglichkeiten der 
Laplace- sowie der Mellin-Transformation auf die Ermittlung von asymptotischen 
Entwicklungen, welche sich meist mit auffallender Leichtigkeit ergeben. Ferner 


; fea nee zs 7 a TA ee fii in ote a 4 
¥s : We ay f ip | me 4 4 
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integralen der Form 


(2) “ar (2) (= 


c 
~ gemacht. . ; kine 
Im zweiten Teil finden wir Anwendungen auf Fakultatenreihen: Auf Grund . 
der Korrespondenz ; Ne 
= = {(1—e*)"} 


s(s + 1) (s+ 2)-+- (s +n) 


- lassen sich verschiedene Aussagen, insbesondere tiber die Ausfiithrung elementarer 
Operationen an Fakultatenreihen, sehr leicht herleiten und die Entwicklung von ~ 
Funktionen nach Fakultatenreihen wird stark erleichtert. 

Im dritten Teil kommen alsdann die bekannten Anwendungen der Laplace- — 
Transformation auf gewodhnliche Differentialgleichungen und insbesondere auf 

_ Probleme der Elektrotechnik zur Behandlung, wobei auch Regelungstechnik und 

Filtertheorie gebiihrend beriicksichtigt werden. Der Verfasser beniitzt hier die in 


es der Elektrotechnik tibliche Terminologie, was dem technisch orientierten Leser 
 — das Studium erheblich erleichtern diirfte. Der dritte Teil befasst sich ausserdem 
- auch mit linearen Differentialgleichungen (und Systemen von solchen), deren 


Koeffizienten variabel, aber analytisch sind. 
Der vierte Teil bringt Anwendungen auf partielle Differentialgleichungen. 
Bi Mit der Laplace-Transformation lassen sich in der Regel nur hyperbolische und 
a; parabolische Typen behandeln, und zwar, indem man die Zeitvariable wegtrans- 
eae formiert. Dabei fiihren eindimensionale Wellen- und Warmeleitungsprobleme 
=e bekanntlich unter gewissen Bedingungen auf gewdhnliche Differentialgleichun- 
a gen, so dass in diesen Fallen eine vollstandige Diskussion der Lésung méglich ist. — 
Der fiinfte Teil ist den Differenzengleichungen gewidmet. Hier ist die Laplace- — 
Transformation dank des Verschiebungssatzes ein besonders bequemes Hilfs- 
ar mittel; freilich nur solange es sich um lineare Differenzengleichungen mit kon-— 
ore, stanten Koeffizienten handelt. Dafiir lassen sich in der gleichen Weise auch Dif- 
; ferentialdifferenzengleichungen behandeln. ; 
cy ae Ebenso ist die Laplace-Transformation fiir die im sechsten Teil behandelten 
iy. Faltungsintegralgleichungen pradestiniert, entspricht doch dem Faltungsintegral 
im Originalraum eine gew6hnliche Multiplikation im Bildraum, und dem komple- 
ee xen Faltungsintegral 
a: 1 
PM 3 


eh he 4 ORE, Se meen 


ft (s — 6) g(a) do 


7 im Bildraum (gleicher Integrationsweg wie beim Umkehrintegral) das Produkt 
F(t) G(¢) im Originalraum. Eine gewisse Schwierigkeit entsteht lediglich dadurch, 
i dass eine Laplace-Transformierte fiir S > oo notwendigerweise gegen 0 konver- 
: gieren muss. 


Schliesslich behandelt der siebente Teil die endliche Laplace-Transformation 
i. b 
my fe 


/ e-st F(t) dt, 


a 


die mit den Funktionen vom Exponentialtypus zusammenhdanegt. 
Bi H. RuTISHAUSER 


zur routinemassigen 
Aufzeichnung von 
Kerninduktions- 
Spektren mit hochster 


Spektrograph a neeaee 


TRUB, TAUBER: ZURICH 
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Einfihrung in Theorie und Anwendung 
der Laplace-Transformation 


Ein Lehrbuch fiir Studierende der Mathematik, Physik und 
Ingenieurwissenschaft. 


Von Prof. Dr. Gustav Dortscu, Professor an der Universitat Freiburg i. B. 


Mathematische Reihe, Band 24 —- Sammlung «Lehrbiicher und Monographien 
aus dem Gebiete der exakten Wissenschaften» 


(1958). 304 Seiten mit 40 Figuren. Ganzleinen Fr. 39.40 (DM 39.40). 


Inhalt (gekiirzt): Das Laplace-Integral von physikalischen und mathema- 
tischen Gesichtspunkten aus — Die Konvergenzhalbebene und die Holomor- 
phie der Laplace-Transformierten ~ Die Abbildung der fundamentalen Opera- 
tionen — Gewohnliche Differentialgleichungen und Systeme von solchen — 
Differenzengleichungen - Das Verhalten der Laplace-Transformierten im 
Unendlichen — Das komplexe Umkehrintegral, Deformation des Weges, Aus- 
wertung des Integrals - Bestimmung der Originalfunktion durch Reihenent- 
wicklung — Die Parsevalsche Gleichung—Asymptotische Entwicklung der Bild- 
und der Originalfunktion — Gewdhnliche Differentialgleichungen mit Poly- 
nomkoeffizienten — Partielle Differentialgleichungen — Integralgleichungen. 


Das dreibandige «Handbuch der Laplace-Transformation» des Verfassers 
stellt das gesamte heute vorliegende Material iiber die Laplace-Transforma- 
tion zusammen und ist in erster Linie als Grundlage fiir die wissenschaftliche 
Arbeit auf diesem Gebiet anzusehen. Im Gegensatz dazu bringt das vorlie- 
gende kurze Lehrbuch nur den Stoff, der fiir den « Normalverbraucher» un- 
erlasslich ist, wenn er die Literatur, die sich der Laplace-Transformation 
bedient, wirklich verstehen oder die Laplace-Transformation bei seinen 
eigenen Arbeiten benutzen will. 
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Jahrbuch 19 57 det Wissenschafilichen Gesellschaft 
fiir Luftfahrt e.V. 


Herausgegeben von Prof. Dr. H. BLEnkK 


Mit den Vortragen der WGL-DVL-Tagung in Essen vom 9. bis 12. April 1957. 
504 Seiten mit 623 Abbildungen und 23 Tabellen. 1958. Leinen DM 58,— 


Ausserdem sind lieferbar: Jahrbuch 1956 DM.38,— 
Jahrbuch 1955 DM 48,— Jahrbuch 1954 DM 36,— 
Jahrbuch 1953 DM 28, — Jahrbuch 1952 DM 24,— 


Eine wichtige Neuerscheinung 
Ballistik 


von Dr. R. E. KuTTERER, 3., verbesserte und erweiterte Auflage. 
Etwa 300 Seiten mit 176 Abbilaiagen! Leinen ca. DM 38, — (erscheint Dezember 1958) 


Viele der in der Ballistik entwickelten Methoden, zum Beispiel die photographischen und kinemato- 

grtaphischen Verfahren, die Druckmessverfahren, die Kurzzeitmessverfahren sowie die Temperatur- _ 

messverfahren, werden heute in der Industrie in steigendem Masse verwendet. Das Interesse an 

diesem Buch diirfte daher auch bei Ingenieuren und Physikern ausserhalb des ballistischen Inter- 
' essengebietes liegen. 


VV Neues Physikalisches Schrifttum aus dem Vieweg Verlag 


Halbleiter und Phosphore 


Vortrage des Internationalen Kolloquiums «Halbleiter und Phosphore» in Garmisch- 
Partenkirchen. Herausgegeben von Prof. Dr, M. Scu6n und Prof. Dr. H. WELKER. 
_ 688 Seiten mit 391 Abbildungen. 1958, Leinen DM 68,— 


detain ht es ts elke 


Dieser Band vereinigt annahernd 100 Vortrage, die von bedeutenden Fachwissenschaftlern aus q ; 


Grossbritannien, USA, UdSSR, Japan, Frankreich, Niederlande, Schweiz, Tschechoslowakei, Polen 
und Deutschland auf dem Internationalen Kolloquium gehalten wurden. Er enthalt auch die Vor- 
trage, die aus Zeitmangel nicht vorgetragen werden konnten. rie 


Halbleiterprobleme 


Herausgegeben und kommentiert von Prof. Dr. Dr.-Ing, e.h. W. ScHoTTKY 


Band IV: Tagung Heidelberg 1957 
376 Seiten mit 75 Abbildungen. 1958. Leinen DM 46,80 


Ausserdem sind lieferbar: 
Band III DM 36,80, Band II DM 28,80, Band I DM 28,80 
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Ne Verlag Friedr. Vieweg & Sohn, Braunschweig 


